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Abstract 



We study static spherically symmetric solutions of Einstein gravity plus an action 

O ■ 

psj ■ polynomial in the Ricci scalar, R, of arbitrary degree, n, in arbitrary dimension, D. The 
^ ■ global properties of all such solutions are derived by studying the phase space of field 
equations in the equivalent theory of gravity coupled to a scalar field, which is obtained 
by a field redefinition and conformal transformation. The following uniqueness theorem is 
obtained: provided that the coefficient a 2 of the R 2 term in the Lagrangian polynomial 
is positive then the only static spherically symmetric asymptotically flat solution with a 
regular horizon in these models is the Schwarzschild solution. Other branches of solutions 
with regular horizons, which are asymptotically anti-de Sitter, or de Sitter, are also found. 
An exact Schwarzschild-de Sitter type solution is found to exist in the R + aR 2 theory 
if D > 4. If terms of cubic or higher order in R are included in the action, then such 
solutions also exist in four dimensions. The general Schwarzschild-de Sitter type solution 
for arbitrary D and n is given. The fact that the Schwarzschild solution in these models 
does not coincide with the exterior solution of physical bodies such as stars has important 
physical implications which we discuss. As a byproduct, we classify all static spherically 
symmetric solutions of D-dimensional gravity coupled to a scalar field with a potential 
consisting of a finite sum of exponential terms. 
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1. Introduction 

Theories of gravity involving higher powers of the Riemann tensor in the lagrangian 
have been proposed in several different contexts since the first days of general relativity. 
They were first introduced by Weyl in his affine theory, which aimed to unify gravity and 
electromagnetism [1] . Such models have become attractive again in recent years following 
the demonstration that the addition to the Einstein-Hilbert lagrangian of terms quadratic 
in the Ricci tensor leads to a renormalisable theory [2]. Unfortunately a massive spin 2 
"ghost" is present in the linearised spectrum, leading to an instability of the theory and 
a loss of unitarity. It has been suggested that this problem would disappear in a full 
non-perturbative treatment of the model, or even that the ghost states in the perturbative 
expansion could be a gauge artifact [3]. However, a non-perturbative formulation is still a 
distant goal, while the latter possibility seems to have been ruled out [4] . 

The effects of higher derivative gravity have also proven to be useful in cosmology, 
beginning with the early work of Starobinsky [5] and Kerner [6] who introduced higher 
derivative terms with a view to obtaining solutions which avoid the initial singularity. 
Later on it was realised that such models can lead to inflationary expansion driven only by 
gravity [7-9]. Higher derivative models have also been studied in the context of quantum 
cosmology [10]. In particular, it has been argued that the introduction of quadratic terms 
may solve some of the problems due to the non-positiveness of the ordinary Einstein-Hilbert 
action for euclidean quantum cosmology [11]. It is also interesting to note that quadratic 
corrections to the action are obtained from quantum wormhole effects [12]. 

Finally, we should mention that higher order lagrangians arise naturally in higher 
dimensional theories, such as Kaluza-Klein and string models. In the first case they 
are introduced in order to obtain spontaneous compactification from purely gravitational 
higher-dimensional theories, and in this context the ghost-free Gauss-Bonnet actions have 
attracted much interest [13,14]. In the second case, they are obtained as an effective low 
energy action [15]. 

While the cosmology of these models has been largely studied, both in four [5-9,16] 
and higher dimensions [17,18], comparatively little is known about black hole solutions. 
The weak field limit has been studied to some extent for the R + aR 2 theory in four 
dimensions [19,20], but the properties of the full solutions in higher derivative theories 
remain largely unexplored. On dimensional grounds, one would expect that the higher 
derivative terms become dominant in the proximity of the singularity. However, it is still 
possible that horizons exist and that, contrary to the assumptions of refs. [19] and [20], 
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analogues of the usual uniqueness theorems for stationary axisymmetric black holes can 
be derived. Indeed, such conclusions can be immediately drawn in the case of the R + aR 2 
theory in four dimensions as a result of the "no hair" theorem proved by Whitt [21], as we 
show in Appendix A. Thus the question of the nature of the static spherically symmetric 
solutions of more general higher derivative models is very interesting from the point of 
view of general relativity. 

One of the reasons for the lack of attention to the black hole problem is undoubtedly 
the difficulty of solving the higher order differential equations arising in higher derivative 
models. Some progress can be made, however, by using the fact that higher derivative 
theories are equivalent, by redefinition of the metric, to ordinary Einstein gravity coupled 
to a scalar plus a massive spin-2 field. This result was first proven by Higgs [22], and 
later rediscovered by Whitt [21], in the case of the R + aR 2 theory in four dimensions, for 
which only the scalar field is present in the effective theory. More recently the equivalence 
has been extended to the case of actions containing powers of the Ricci and Riemann 
tensor [23] The importance of this equivalence is that one can now use the formalism of 
ordinary general relativity to study the more general higher derivative theories. 

Our analysis in this paper will make use of the equivalence of the general D- dimen- 
sional action 

X An 2 1 J 



to Einstein gravity coupled to a scalar field [23]. Here / is an arbitrary function of the 
Ricci scalar R, and k 2 denotes the gravitational constant in D dimensions. If we define a 
by 

where 



In [4eK 2 f'(R)] , (1.2a) 



e=| 1 ' if//>0 ' (126) 

£ \-i, iff <o, [L - Zb) 

and make a conformal transformation 

g ab =[4eK 2 f'(R)f {D - 2) g ab: (1.3) 
then the field equations derived from (1.1) are equivalent to those derived from the action 

s = I d D x^§ | A _ -L_g^ daadba _ v{a) 1 j {1Aa) 



For a different approach to the problem see [24] . 
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where 



V=— 2 [Ae^f\R)]- D/{D - 2) [Rf'-f] 



An 2 



(1.46) 



with R, f and /' denned implicitly in terms of a via (1.2). For the quadratic theory, for 
example, with / = R + aR 2 we find [18] 



V = 



/ 2(D-4)k<7 
exp ' 



16k 2 



a 



1 — e exp 



-2kct 



Similarly, for the cubic theory with / = R + aR 2 + bR 3 we find [16] 



V = 



54n 2 b 2 



exp 



-2Dkct 



(D-2)y/D=T 



± 



a 2 — 36 ( 1 — e exp 



9 



2kcj 



3/2 



a H — a& 1 — eexp 



2kct 



For other actions polynomial in 



(1.5) 



(1.6) 



f(R) = R + J2%R P > 

p=2 



(1.7) 



and /' can be inverted to give a precise analytic expression for R in terms of a, and hence 
for V(cr), for a general polynomial only if n < 5. 

The general higher order field equations obtained from the action (1.1), with / of the 
form (1.7), are given by 



1 



Rab ~ ~9abR + J2 a P { P RP ~ lR ab ~ P(P ~ 1)# P " 3 RR;ab + (p ~ 2)R;aR;b 



p=2 



+ 9ab 



pip - l)R p ~ 3 (rDR + (p - 2)K c R,}j - X - 



R p 



= 0. 
(1.8) 



Since the arbitrary-dimensional Schwarzschild solution 



ds 2 = - 1 



- 75=3 J dt ^s) dr + dQ n-2, (1-9) 
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has R = throughout the domain of outer communications, it is clear that the Schwarz- 
schild solution solves the equations (1.8) for any choice of the constants a p . Thus the 
problem before us is to determine whether the Schwarzschild solution is the only static 
spherically symmetric asymptotically flat solution with a regular horizon. On account of 
the equivalence of the higher order theory to the theory described by the action (1.4) this 
uniqueness problem can be regarded to as the problem of establishing a "no hair theorem" 
for the latter model. 

Whitt established such a theorem in the case of the R+aR 2 theory in four dimensions^ 2 
[21], by demonstrating that all asymptotically flat stationary axisymmetric solutions to the 
higher order vacuum equations must have R = in the domain of outer communications. 
If one adds extra matter fields to the action one still finds that R = for stationary, 
axisymmetric, asymptotically flat solutions provided that the energy-momentum tensor is 
traceless and satisfies the matter circularity condition. For such solutions, therefore, (1.8) 
becomes equivalent to the usual Einstein equations, and the usual uniqueness theorems and 
no hair theorems will carry over to the fourth order theory. If one considers an arbitrary 
polynomial in R of the form (1.7), however, then Whitt' s argument breaks down, as we 
demonstrate in Appendix A. Consequently a different approach is called for. 

Our approach here will differ not only from that of Whitt, but also from other standard 
approaches to no hair theorems [26,27], in that we will solve the problem by studying the 
phase space of the field equations obtained from (1.4). We will take advantage of the fact 
that by making a judicious choice of coordinates these equations may be written in the 
form of a 5-dimensional autonomous system of ordinary first-order differential equations, 
so that all the global properties of the solutions can be derived. Our approach not only 
has the advantage that it can be used to establish a black hole uniqueness theorem for 
general actions polynomial in R, but it will also enable us to determine the nature of other 
solutions in these models which have regular horizons but which are not asymptotically 
flat. 

The analysis we will use here is very similar to that developed in refs. [28] and [29], 
henceforth denoted I and II respectively, where we derived the global properties of static 
spherically symmetric solutions in models of gravity which arise from the dimensional 
reduction of certain higher-dimensional gravity theories. The scalar field in (1.4) then 
corresponds to the radius of the extra dimensions, (the "compacton" ) , and the potential 

^ 2 In fact, a black hole uniqueness theorem for the more restricted case of pure B? theory, without an 
Einstein-Hilbert term, was obtained much earlier on by Buchdahl [25]. 
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V(cr) contains one or two exponential terms. 

The analysis of I and II is based on the fact that the appropriate field equations 
can be reduced to a 5-dimensional autonomous system of first-order ordinary differential 
equations. This is possible essentially due to the fact that the field equations form a system 
very similar to those of a Toda lattice when written in terms of appropriate coordinates 
[30], (the Toda lattice being an integrable system). Since the metric and scalar fields 
are related to the functions X, Y, V, Z and W of the 5-dimensional phase space, Ai, 
they are necessarily regular at all points of the integrals curves apart from critical points. 
Consequently, in order to determine the global properties of all solutions - namely the 
structure of their singularities, horizons and asymptotic regions - it suffices to study the 
properties of the solutions at critical points of A4. In order to determine which critical 
points are connected to which other ones by integral curves requires a careful analysis of the 
structure of the space M.: in particular of surfaces corresponding to particular subspaces, 
which separate integral curves corresponding to spacetimes of different causal structures. 
This method represents a powerful analytical tool for these (and possibly other) models 
for which it is not possible to write down the general static spherically symmetric solutions 
in a closed analytic form. 

One should mention that on account of the way in which the phase space functions 
are constructed, the metric functions (with signature — (- + •••+) are necessarily positive. 
Thus for the Schwarzschild solution the integral curves obtained correspond to the domain 
of outer communications only. A similar analysis of the Reissner- Nordstrom solution yields 
a phase space with distinct regions which correspond to (i) integral curves in the domain 
of outer communications; and (ii) integral curves in the region between the singularity and 
the inner horizon. (These two regions of the phase space are separated by a surface which 
corresponds to the Robinson-Bertotti solutions.) Thus the method described does not pick 
out regions of the spacetimes in which the Killing vector d/dt is spacelike. Such regions 
can only be obtained by continuation of such solutions as are described here. 

The first step of our analysis here will be to generalise the work of I and II to establish 
the properties of solutions to the equations derived from (1.4) with a potential consisting 
of a sum of s exponential terms 

V(,) = ^ (MO) 

1=1 X 7 

the Aj and g v i = l,...,s being constants. Of the higher derivative theories, only the 
R + aR 2 theory has an action which is precisely of this form. However, for each n there 
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will be a special choice of the constants a p for which the potential (1.4b) reduces to this 
form. This will be so if the (n — l)-th root of /' is linear in R. In particular, if we choose 

(2a 9 )P- 1 (n-2)\ 

°* = (n-l)P-* P !(n- P )! ' 3 " P ~ (LH) 



then 



n-l 



and we find that the potential takes the form 



(n-l) 2 / —2Dna \ I / 2nna 
V = „ — exp ^=^= < e exp 



8nn 2 a 2 \(D- 2)y/D -1/1" \(n - - 1 



exp 



where^ 3 



n-l 



1, if 1 + 2a 2 R/(n - 1) > 
-1, if 1 + 2a 2 R/(n - 1) < 




We shall call this class of models "special polynomial R theories" . 

Before proceeding further, we note that the definition of asymptotic flatness in our 
analysis requires some clarification. On account of the conformal transformation (1.3) the 
radial coordinate of the higher derivative theory, r, is related to the radial coordinate of 
the effective theory, f by 

r = [toV(fl)] - 1/( °- 2) f = exp ( - _-^ o = ) r . (1.14) 

Thus solutions which are asymptotically flat in the higher derivative theory need not be 
asymptotically flat in the effective theory, since it is the dependence on r as measured 
by the metric g a b which is of physical importance rather than the dependence on f as 
measured by the metric g a b- Moreover, the definition of the asymptotic region can vary 
between the physical and the effective theories. The definitions of spatial infinity and of 



^ 3 Note that e = e n ~ l with these definitions: for n even e and s have the same sign. However, if n is 
odd then s = +1, and the sign of s corresponds instead, for example, to the plus or minus sign multiplying 
the first term in (1.6). 
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asymptotic flatness will only coincide when a — * const, as f — > oo. However, since a is not 
a physical field here, (unlike the Kaluza-Klein case [28,29]), we need place no requirements 
on its asymptotic form at spatial infinity in the effective theory. Thus a discussion of 
the existence of static spherically symmetric black holes in the higher derivative models 
requires an examination of all static spherically symmetric solutions in the effective theory, 
and not merely the ones which are asymptotically flat according to the effective theory. 

We will begin by deriving the global properties for the theory with the potential (1.10) 
- almost all cases are dealt with in §2. For special values of the parameters, which include 
the cases of the R + aR 2 theory and all potentials (1.13), the structure of the phase space 
is modified. Such solutions will be discussed in §3. We will then generalise the analysis to 
include the potential (1.6) to that of the R + aR 2 + bR s theory in §4. This generalisation 
will be found to produce only minor changes to the results of §2 and §3. We will further 
demonstrate that our results can be extended to the case of the polynomial of arbitrary 
degree. Some physical implications of our work are discussed in §5. 



2. The general exponential sum potential 



2.1 The dynamical system 

As in I and II we will begin by choosing coordinates 

g ab dx a dx b = e 2il (-dt 2 + f 2m d£ 2 ) + f 2 g a pdx a dx fi , (2.1a) 

where m = D — 2, u = u(£)> ^ = ^(0 an d 9ab is the metric on an arbitrary m-dimensional 
Einstein space: 

Rap = (m - l)\g a p- (2.16) 

Of course we are most interested in the case in which g ab is the metric on a 2-sphere (and 
A = 1). However, the structure of the phase space is better revealed by taking the more 
general ansatz (2.1b). In order to write the equations as a first order system which is 
everywhere well-defined, we will choose an ordering of the g. such that 



gl >g 3 >g 4 > ... >g s >g 2 . 



(2.2) 
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If we now define the functions" 4 £, rj and \ by 

( = u+ (m- l)lnf, (2.3) 

?7 = w + mlnf-^ — , 2.4 

m 

X = u + m\nf . (2-5) 

m 

then the field equations become 

C" = (m - l) 2 Ae 2 ^ + A^ 2 " + A 2 e 2X + ^ V 2 ^ X+ ^), (2.6a) 

77" = rn(jn - 1) Ae 2C + — ( m + 1 - g A \ e 2r] + — (m + 1 - g, g 9 ) A. 

/// V 1 / 1 /// ~ 

+ l^( m+ l-g lg J A ,e 2 (- X +^), 

i=3 

x " = m {m - l)Ae 2C + — (m + 1 - gl g 2 ) A^ 2 " + — ( m + 1 - g 2 ) A. 
with the constraint 



. e 2X 



(2.66) 



. e 2X 



(2.6c) 



m 

i=3 



, , lV /2 , 2mC(g 2 ^-g lX / ) , (l + (m-l)g 2 2 ) /2 (l + ( m -l) gl g 2 ) y , 

(m + l)C + + \2~^ ~ 2 7 \2~^ ^ 

gi~g 2 (gi-g 2 ) (gi-g 2 ) 

(l + (m - l) gl 2 )^ /2 _ ^ c + ^ + A, g2X + J_ ' Xie2(aiX+ p iv) = Q 

™ m m 



(g 1 -g 2 ) 2 " " " m 



i=3 



(2.6d) 

where, for s > 3, 

a . = g i ^g i) i = 3,...,s, (2.6e) 

Si §2 



and 



# = l- ai = -?i— §2, i = 3 ,... jSj (2.6/) 
Si-g 2 



#4 



These functions correspond, in fact, to the differences of the Toda lattice coordinates [30]. 
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If s < 3 then the last summation term in (2.6d) vanishes. Note that on account of (2.2), 
< ai < 1, < fa < 1, and aij = Pi only in the special case in which s = 3 and 

g 3 = (Si + S 2 ) /2- 

These equations can be recast in the form of a 5-dimensional autonomous system of 
first-order differential equations. If we define variables V, W, X, Y and Z by 

V = X ', W = e*, X = C, Y = V ', Z = e\ (2.7) 

then the constraint (2.6d) can be regarded as a definition of e 2C - . Eliminating the e 2C - terms 
from (2.6a-c) we therefore obtain the system " 5 

X' = 1 ( A 1 Z 2 + X 2 W 2 + y XiW 2ai Z 2 ^ ] - (m ~ 1)P (2.8a) 

m I U J 

Y' = ^ | (g/ - l) \Z 2 + ( gl g 2 - 1) X 2 W 2 + £ (g lg . - 1) \W 2 ^Z 2 ^ | - P, (2.86) 

^ = ^ | (SiS 2 - 1) \^ + (g/ - l) A 2 ^ 2 + £ (g 2 g, " 1) \W 2a% Z 2<3% | - P, (2.8c) 

Z' = YZ, (2.8d) 
W = VW, (2.8e) 

where 



P = m<( (m + 1)X 2 + ••— 2L_^ + 



2m^ 



:(g 2 y- gl F) (i + (m-i)g 2 2 ) 



gi - s 2 (gi - g 2 ) 2 

(2-8/) 

2 (1 + (m - l) gl g 2 ) (l + (m - l) gl 2 ) 

(gi-g 2 ) 2 (gi-g 2 ) 2 

The dynamical system of the Kaluza-Klein models dealt with in II, with n e extra 
dimensions and a higher-dimensional cosmological constant A, is retrieved by setting s = 2 
and 

gi = g2 -i = v^_+^e > A 1 = n e (n e - 1)A, A 2 = -2A. (2.9) 



»5 



We are using a different normalisation here for Z and W compared with that used in I and II. 



- 11 - 



In the case of the R + aR 2 theory, and indeed any higher order theory given by the 
potential (1.13), the appropriate dynamical system is obtained by taking equations (2.8) 
and setting s = 3, 

m + 2 2(n-l)-m 1 , 

g i = o /—tt > g 2 = W< TV^^f ' g s = r~^T > (2 - 10a) 

2\/m + 1 2{n — l)y/m+l v m + 1 

n—1 1 
a 3 = — , Pa = (2.106) 



n n 



and 



1 2na 2 2 2na 2 3 2a 2 v ' 

As in the case of the Kaluza-Klein models the phase space has a great many symmetries 
which greatly simplify the analysis. Equations (2.8b) and (2.8e) ensure that trajectories 
cannot cross either the W = or the Z = subspaces. These two subspaces correspond 
physically to the cases in which A 2 = and A 1 = respectively. As we have written them, 
equations (2.8) are valid for W > and Z > 0. It is possible to make the equations valid 
for all W and for all Z by introducing modulus signs in the terms W 2cti and Z 2 @ i which 
involve fractional powers of Z and W. However, this merely introduces a trivial symmetry 
between trajectories in the W > and W < portions of the phase space, and between 
trajectories in the Z < and Z > portions of the phase space. Thus we may restrict our 
attention to Z > and W > without loss of generality. 

The hyperboloid defined A = 0, or 

s 

P + \Z 2 + X 2 W 2 + XtW 2 ^Z 2 ^ = 0, (2.11) 

i=3 

similarly forms a surface which trajectories cannot cross. It partitions the phase space into 
the two physically distinct regions with A > and A < 0. 

If W = and g 1 ^ then 



_ gi(gi ~ § 2 )( toX + c i) + (1 + (m - l)g 1 g 2 ) Y 
1 + (m - l)g^ 



v = 6iV&i 62A ; ; ^ \^ 2 V J&i&2J , (2.i2) 



while if Z = and g 2 ^ then 



v ~g 2 (gi - S 2 )(mX + c 2 ) + (1 + (m - l) gl g 2 ) V 

1 + (m - l)g 7 ' (2 - 13) 
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where c ± and c 2 are arbitrary constants. (If g 1 = or g 2 = we have instead V = 
Y + const.) Physically (2.12) is equivalent to s = 1, (i.e., A^ = 0, i > 2), while (2.13) gives 
the same system with A 1 — > A 2 and g 1 — > g 2 . Thus in each case a further degree of freedom 
can be integrated out, giving rise to a 3-dimensional autonomous system. The properties 
of such systems were studied in I and II. Further simplifications arise if in addition one of 
the constants A or A 1 (or A 2 as appropriate) is zero. In these cases it is in fact possible to 
integrate the field equations exactly. For completeness we list these solutions in Appendix 
B. 



2.2 The (anti)-de Sitter subspaces 

In addition to the W = 0, Z = and A = subspaces, there is at least one other 
3-dimensional subspace which exists for all s > 2, which was not studied in detail in II. 
These subspaces may be identified by noting that solutions for which a is constant globally 
form a special class. The actual value which this constant takes may be determined from 
the field equations. In particular, we find that if V = Y and W = 7Z, (7 > 0), where 

s 

\Si + A 2 S 2 7 2 + E X ^il 2ai = 0, (2-14) 
then the field equations (2.8) reduce to the 3-dimensional system 

X' = -AZ* - (m ~ 1)P (2.15a) 

m 

Y' = -AZ 2 - P, (2.156) 

Z' = YZ, (2.15c) 

where 

A=^(A 1 + A 2 7 2 + E^7 2ai V ( 2 - 15rf ) 



i=3 



and P now simplifies down to 

P = m(X -Y)((m + l)X - (m-l)Y) . (2.15e) 
In the case s = 2, for example, (2.14) is satisfied by any g 1 and g 2 provided that 
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is real. For s > 2, (2.14) may have more than one solution. 

If s = 1 we of course obtain equations (2.15) if g x = and A 1 = -mA. Consequently, 
equations (2.15) are formally equivalent to those of the Schwarzschild-de Sitter solution 
in (m + 2) dimensions, and can be integrated explicitly in terms of the coordinate f. We 
therefore obtain the solution 

g ab dx a dx b = -Adt 2 + A^df 2 + r 2 g a/3 dx a dx^ ', (2.17a) 



with 



C A^- 2 Si/(gi-g2)f 2 

A = A —— - — -, (2.176) 

fm-1 m+1 



and constant scalar field 

e 2„a = 7 m/(g 1 -g 2 ) j (2.17c) 

C being an arbitrary constant. In the case of the higher order theories with potentials 
(1.13), a conformal transformation back to the original metric yields the solution 

ds 2 = -Adt 2 + A _1 dr 2 + r 2 g af3 dx a dx p , (2.18a) 



with 



A = A - - ^ , (2.186) 

r m-l m+1 



and M = C y[2G 7 2 ("»- 1 )(»- 1 )/("»»)]. (We have rescaled t in obtaining (2.18).) The solutions 
are asymptotically de Sitter if A > (or anti-de Sitter if A < 0). In the case of the R + aR 2 
theory, we have 

- 2GM mP 

A - A - + (m + 1)(m + 2)(m _ 2)a , (2-19) 

and consequently solutions of this type exist only for m > 2, i.e., for D > 4. These 
solutions are asymptotically anti-de Sitter if a > 0. For n > 2 solutions of the type (2.18) 
also exist for D — 4. In general, it is possible to find more than one branch of solutions - 
some of which are asymptotically de Sitter, and some of which are asymptotically anti-de 
Sitter. 

If M = we retrieve cosmological de Sitter solutions, the existence of which has 
been discussed previously by Barrow and Ottewill [8] for an arbitrary f(R) lagrangian in 
four dimensions, and by Madsen and Barrow [31] for more generalised higher derivative 
lagrangains in arbitrary dimensions. 
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The solutions (2.17) do not exhaust all possible solutions to equations (2.15). In par- 
ticular, a special class of solutions for which r(£) is everywhere constant are also admitted. 
These solutions may be determined by direct integration of (2.15) in the case that Y = X. 
Since P = in this instance we have Y' = X' and (2.15) reduces to a 2-dimensional au- 
tonomous system. The equations can be readily integrated using the coordinate Z = e v , 
and we find 

g ab dx a dx b = y*i/fei-*2) |-Z 2 dt 2 + + (^^) gaffiUc"*?} , (2-20) 

where C is an arbitrary constant, and Z > 0. (We have used the freedom of rescaling 
t to remove an unphysical constant.) Thus we find a solution with A > provided that 
A > 0. These solutions are topologically a product of 2-dimensional anti-de Sitter space 
with an m-sphere of constant curvature - a type of Robinson-Bertotti solution. Solutions 
with A < and A < are similarly a product of 2-dimensional de Sitter space with 
m-dimensional hyperbolic space. 

To compare these results with our later studies of the 5-dimensional phase space it is 
useful also to give a description of the 3-dimensional phase space here. 

Provided that A^O then the only critical points at a finite distance from the origin 
are the lines 

Z = 0, Y = X, (2.21) 



and 



Z = 0, F= (m+1)X (2.22) 
m — 1 



These solutions have A = as well as W = Z = 0. Solutions lying in the Z = plane, as 
depicted in Fig. 1, are Schwarzschild solutions. Equations (2.15) can be integrated directly 
in this case since 

777 

Y = -(X + k), (2.23) 

m — 1 

where k is an arbitrary constant. The critical points in the first quadrant are found to 
correspond to the limit £ — > — oo, while the critical points in the third quadrant correspond 
to the limit £ — > oo. We find that endpoints of trajectories on the line Y = X correspond to 
horizons with f — > const., while endpoints of trajectories on the line Y = (m+l)X/ (m — 1) 
correspond to r — > singularities. The Z = trajectories ending on the line Y = X 
thus represent positive mass Schwarzschild solutions, while those ending on the line Y = 
(m+l)X/(m — l) represent negative mass Schwarzschild solutions with naked singularities. 
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The constant k in (2.23) is related to the Schwarzschild mass: the k = trajectories 
represents flat space. Furthermore, if a trajectory which corresponds to a positive mass 
Schwarzschild solution in one quadrant is traced back to the opposite quadrant, then one 
obtains the negative mass Schwarzschild solution of the same absolute mass in the A > 
region. 

If A = then additional critical points exist near the origin. Furthermore, (2.15) 
reduces to a 2-dimensional autonomous system for all X, Y and Z. We will defer discussion 
of this special subspace until §3. 

Small perturbations about the critical points (2.21) and (2.22) yield the eigenvalues 
A = 0, Yq, 2Xq in the 3-dimensional subspace, (where Yq = Xq or Yq = (m+l)Xo/(m— 1) as 
appropriate), and additional eigenvalues A = 0, Y Q for perturbations in the extra directions 
in the full 5-dimensional phase space. The pattern of trajectories is therefore identical 
to that of the corresponding trajectories in I and II. In the 3-dimensional subspace each 
critical point (X o ,Yq,0) in the first (third) quadrant repels (attracts) a 2-dimensional set 
of trajectories which lie approximately in the plane Y = (ml =F Xo)/(m — 1). Since the 
trajectories are approximately planar the one zero eigenvalue corresponds to the degenerate 
direction perpendicular to this plane. In the 5-dimensional phase space there is an extra 
degenerate direction, and the dimension of the set of trajectories repelled (attracted) for 
first (third) quadrant trajectories increases by one. 

To complete the description of the phase space it is necessary to describe the critical 
points at infinity. Following I and II, we introduce spherical polar coordinates 

X = p sin 9 cos <j> Y = p sin 9 sin Z = pcos9. (2.24) 

The surface at infinity is then brought to a finite distance from the origin by the transfor- 
mation 

p = p(l-p)-\ 0<p<l. (2.25) 

If we define a coordinate r by dr = pd£ = p(l — p) _1 <i£, then on the sphere at infinity, i.e. 
at p = 1, dp/dr = identically while 



d9 

— = cos 9 1 —A cos 2 9 [cos <fr + sin <p] — sin 2 9 
dr 1 



. , = f(m- 1) , . A 
sin 6 + P. [ cos 6 + sin 6 

V m J 



(2.26a) 

d(j) 1 



dr sin 9 



-A cos 9 (cos 4> — sin 4>) + sin 9P 1 



(m - 1) . 

sin <p — cos ( 

TO 



(2.266) 
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where 

P 1 = m ((m + 1) cos 2 (p — 2m cos sin + (m — 1) sin 2 0) . (2.26c) 
Four sets of critical points are found: 

(%). First of all we obtain the endpoints of the lines of critical points Y = X and Y = 
(m + l)X/ (m — 1), for which A = W = Z = 0. The points located at 

e=l (2.27a) 

or X = ±oo, Y = X, Z = 0, (2.276) 

will be denoted L\ and L 2 . The points located at 

^=2' ^ = arctan (^ — J J (2.28a) 

or X = ±oo, Y = (j^-j^J x , Z = Q - ( 2 -286) 

will be denoted L3 and L4. As for points at a finite distance from the origin, the points 
Li and L2 correspond to horizons, and L3 and L4 to f — > singularities. The points L13 
(-^2,4) repel (attract) a 2-dimensional set of trajectories, which are unphysical, however, 
since they are confined to the sphere at infinity. 

(ii). Two critical points, which we will denote Mi and M 2 , are located at 

7T / in 

6= -, (j> = arctan ) (2.29a) 

2 \m — 1 

(vn \ 
)X, Z = 0, (2.296) 
m — 1 J 

in the A > portion of the phase space. These points correspond to the asymptotic region 
(f — > 00) of the Schwarzschild solutions which lie in the Z = plane. These are the only 
trajectories which end on these points: they are found to be saddle points with respect to 
other directions in the phase space. 
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(iii) . If A < then there are two critical points, which we will denote S\ and S2, which 
are located at 

6 = arctan V-2A, (p = - , — (2.30a) 

or X = ±00, Y = X, Z = —=, (2.306) 

V-A 

in the A < portion of the phase space. At these points we find that f —> const. These 
points correspond to the Z — > 00 region of the Robinson-Bertotti-like solutions (2.20). 
An analysis of small perturbations in the 3-dimensional subspace reveals that the point 
Si (S2) in the first (third) quadrant is an attractor (repellor) for a 2-dimensional set of 
trajectories: these are of course the solutions (2.20) with A < and A < 0, which lie in 
the plane Y = X. This clarifies the nature of the corresponding points Sis m H, which 
were not discussed in detail there* 16 . 

If A > (and A > 0) then no Z — > 00 region is defined since the maximum value Z 
can take is (C/A) 1 / 2 . Thus instead of ending on the points Si^ trajectories move from 
one point on the Z = 0, Y = X line to another point on the same line in the opposite 
quadrant. This pattern is made clear by Fig. 2 where we plot the Y = X plane through 
the 3-dimensional subspace. 

(iv) . If A < then there are two critical points, which we will denote T\ and T 2 , which 
are located at 

/-(2m 2 + 2m + 1)A\ 1/2 , fm + l\ , rt oi . 

= arctan , 6 = arctan 2.31a 

V m+ 1 J \ m J ' 

X = ± oo, y = <"' + 1 >* Z = ±(^±±Y\ (2.316) 

m m \ —A / 

in the portion of the phase space with A = 0. Point T\ (T 2 ) in the first (third) quadrant is 
found to attract (repel) a 3-dimensional set of trajectories in the 3-dimensional subspace. 
These points of course correspond to the asymptotic (f — > 00) region of the Schwarzschild- 
anti-de Sitter solutions (2.17) if A > 0. T 12 are endpoints for trajectories both in the A > 
and A < portions of the phase space. 

If A > then no asymptotic region is defined. Instead the trajectories move from 
the Y = X or Y = (m + l)X/ (m — 1) line in one quadrant to one of these two lines in 



^ 6 Since we are restricting the analysis to W > and Z > here, we of course obtain half or quarter 
as many critical points as in I and II for those points with IV / or Z / 0. 
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the opposite quadrant. Trajectories starting and finishing on the Y = X line with A > 
correspond to the region of the positive mass Schwarzschild-de Sitter solution between the 
Schwarzschild and de Sitter horizons. Trajectories with endpoints on both the Y = X and 

Y = (to + l)X/(m — 1) lines represent the region of the negative mass Schwarzschild-de 
Sitter solution between f = and the de Sitter horizon. Trajectories cannot have two 
endpoints on the Y = (to + l)X/ (to — 1) line since they cannot cross the Y = X plane. 

The pattern of trajectories on the hemisphere at infinity of the 3-dimensional subspace 
is sketched in Fig. 3. Although these trajectories are unphysical it is helpful to sketch them 
since by continuity arguments they will determine the behaviour of the physical integral 
curves which lie within the hemisphere at infinity but near its surface. 

2.3 The W = and Z = subspaces 

These subspaces were discussed in I and II for particular values of g 1 and g 2 . The 
properties for general g ± and g 2 are similar. If A = the only critical points at a finite 
distance from the origin in the full 5-dimensional phase space have both W = and Z = 0, 
and so are common to both subspaces. The critical points are located at X = X and 

Y = Yq. where 

l*ol > ( tTT—\ 2 ) ISi c il' ( 2 - 32 ) 
yl + (m- l)gi y 

and Yq is given by solving the quadratic equation 

( m _ i)y fl 2 _ 2mX Y + (to + 1 - gl 2 )A 2 + gfc? = 0, (2.33) 

with V given by (2.12) with X = X . 

Consider the W = subspace. If g 2 < to + 1 then the pattern of the trajectories is 
the same as in I and II since all critical points lie in the first and third quadrants. Points 
with Y > (first quadrant) correspond to the limit £ — > —00, and those with Y < 
(third quadrant) to the limit £ — > +00. Each point in the first (third) quadrant repels 
(attracts) a 2-dimensional set of trajectories which lie approximately in the plane Y = 

1/2" 



(^l) [ X ± ((! + (rn - l) gl 2 ) A 2 - (to - 1)^) 



There is a zero eigenvalue 

corresponding to the degenerate direction perpendicular to this plane. In terms of the 
coordinate f one finds that f — * at all critical points except those for which c ± = mk, 
which correspond to horizons. These special critical points are of course those lying in the 
W = 0,Z = 0,V = Y plane (c.f. §2.2). Fig. 1 thus represents the plane which bisects the 
W = 0, Z = subspace to pick out the Schwarzschild solutions. 



If g 2 > m + 1 then points for which \X \ < X 1: where 




\Xq\ > Xi, on the other hand, then one critical point lies in each quadrant. The two points 
in the first and third quadrants have the same properties as before. At the critical points 
in the second and fourth quadrants we find that r — > oo. As before there is one degenerate 
direction corresponding to a zero eigenvalue. However, both points are now saddle points 
with respect to the remaining two directions. Each critical point in the second (fourth) 
quadrant repels (attracts) one trajectory from the A > region and one trajectory from the 
A < region of the phase space: these are the Z = solutions discussed in Appendix B.l, 
and depicted in Fig. 4(b). Each point in the second (fourth) similarly attracts (repels) one 
trajectory for each sign of X x : these trajectories are in fact the A = solutions discussed 
in Appendix B.2, and depicted in Fig. 5(c). Their asymptotic form is given by (B.27). 
Thus trajectories in the A > and A < regions of the phase space for which Z is not 
identically zero do not have endpoints in the second and fourth quadrants. If such solutions 
have asymptotic regions, the limit f — > oo must be approached at critical points at the 
phase space infinity. 

The phase space infinity of the W = subspace can be studied once again by intro- 
ducing coordinates (2.24), (2.25). The following equations are obtained for the angular 
coordinates on the p = 1 sphere at infinity 




(2.35a) 




ch/> 
d7 



sin I m 




— cos 



2 9 ((1 -g^cos^-sin^J +sin 2 #P 2 



(m — 1) 



rn 



sin 4> — cos <f> 




where 




1 + (m - l)g : 2 



rn 




m + 1 - g 



cos 2 (j) — 2m cos <fi sin <fi + (m — 1) sin 2 <fi \ . (2.35c) 



Four sets of critical points are found: 
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(%). Once again we obtain the endpoints of the curves of critical points with A = W 
Z = 0. These points, which we will denote L 5 _ 8 are located at 



4> = arctan 



m± y /l + g 1 2 
to — 1 



(2.36a) 



or 



X = ±oo, 



Y = X 



m± y / l±g 1 2 



to — 1 



Z = 0. 



(2.366) 



These points have the same properties as those outlined above for the appropriate cases 
of the A = W = Z = points at finite distances from the origin. 

(ii). We once again obtain the critical points M\ and Mi given by (2.29). As before, these 
points, corresponding to asymptotically flat solutions, are found to act as saddles with 
respect to all trajectories other than the 2-dimensional bunch ofW = Z = Schwarzschild 
solutions. 

(Hi). If A 1 > and g 2 < m ± 1 or if A 1 < and g 2 > m ± 1 then there are two critical 
points, which we will denote N\ and N 2 , which are located at 



6 = arctan 



A.| m+1 " g ' 2 + '" 



-,1/2 



m 



m + 1 - g x 2 



= arctan 



m ± 1 - g x 2 



m 



or 



X = ±oo, Y 



o \ 1/2 

m±l — g/\ fm + l — g, 



m 



mA, 



on the A = surface. In the full 5-dimensional phase space these points also have 



W = 0, 



V = 



(to + 1 - g x g 2 ) 



X. 



TO 



(2.37c) 



If gj 2 < 1 then the point N\ (AT 2 ) in the first (third) quadrant attracts (repels) a 3- 
dimensional set of trajectories in the 3-dimensional subspace. (Ni attracts all trajectories 
in the A > region, and some A < trajectories if g^ < 1.) If 1 < g^ < to ± 1 then N± 
(N 2 ) only attracts (repels) the 2-dimensional set of trajectories lying in the A = surface, 
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and acts as a saddle with respect to other directions. If g x 2 = m + 1 then Ni and -/V 2 
are degenerate with the points L5 and L7. If g^ 2 > m + 1 then the point N\ (N 2 ) lying 
in the fourth (second) quadrant repels (attracts) a 3-dimensional set of trajectories in the 
subspace. 

(iv). If A 1 > then there are two critical points, which we will denote Pi and P 2 , which 
are located at 

1/2 



= Cretan (l + (m - l^jj , = -, — 



(2.38a) 



or X = ±00, Y = X, Z= I ? — r I X (2.386) 



\ (l + (m-l) gl 
In the full 5-dimensional phase space these points also have 




1 + (m - l)g^ 



W = 0, F = ( - 7 ^ &1 7^ 2 ] X. (2.38c) 



If g 2 < 1 these points lie in the A < portion of the phase space, while if g > 1 they 
have A > 0. If g 2 = 1 they are degenerate with points Ni and N 2 . 

If g 2 < 1 the point Pi (P 2 ) in the first (third) quadrant attracts (repels) a 2- 
dimensional set of A < trajectories, acting as a saddle with respect to other directions. 
The 2-dimensional separatrix separates A < trajectories with an endpoint on Ni (N 2 ) 
from trajectories with two endpoints on the A = W = Z = curve. If g^ > 1 then Pi 
(P 2 ) attracts (repels) a 3-dimensional set of trajectories: all A > solutions apart from 
those lying in the Z = plane. 

In Figs. 6 and 7 we sketch trajectories on the hemisphere at infinity of the W = 
subspace for the various cases which give distinct behaviours. 

The properties of the Z = subspace follow by symmetry upon making the substitu- 
tions Y <-> V, Z <-> W, X 1 — > A 2 and g x <-> g 2 in the above discussion. To set our notation, 
the critical points on the sphere at infinity located at 



X = ±oo, V = 



z = 




(2.39) 



- 22 - 



will be denoted R\ and R 2 . The points located at 

1/2 

X = ±00, V = X, W = I t — r- I X, 



A 2 (l + ( m -l)g 2 2) 



(2.40) 



^ = 0, y= I t^7 1* 

l l+( m -l) g 2 



will be denoted Qi and Q 2 - 



2.4 Global properties of solutions 

We turn now to the global properties of the solutions as deduced from the nature 
of trajectories in the full 5-dimensional phase space with A = 0. The behaviour of the 
trajectories may be pieced together in a relatively straightforward manner from the prop- 
erties of trajectories in the subspaces already discussed, since if A = the only critical 
points at infinity other than those already found are the extension of points Lis to the 
one-parameter family of critical points which coincide with the intersection of the A = 0, 
Z = 0, W = surface and the sphere at infinity. We shall denote the whole set {L(y)} : 
where 



m 



- + (m - l) gl 2 ^ m + + (m - 



< y < i • (2.4i) 

m — 1 m — 1 

The points L(y) are located at 

X = ±oo, Y = yX, Z = 0, W = 0, 

( mg 1 {g 1 - g 2 ) + (1 + (m - l)g 1 g 2 ) y ± yj2my - (m - l)y 2 + - m - l\ 
V = [ 1 + (m - l) gl 2 J X - 

(2.42) 

Generically, apart from a few exceptions^ 7 trajectories which are not confined to the 
sphere at infinity have at least one critical point on the A = W = Z = curve. Provided 
that g l 2 < m + 1 and g 2 2 < m + 1, then all such critical points take values X , Y and 



»7 



The W = 0, A = 0, g 2 > 1 solutions given by equations (B.29)-(B.35) in Appendix B have endpoints 



on A?i and N2 on the sphere at infinity, and have no asymptotic region. A similar class of solutions exists 
for Z = 0. Also, in the case of the Robinson-Bertotti type solutions (2.20), some trajectories join the 
points Si and S2 (c.f. Fig. 2(b)). 
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Vo which are either all positive, or all negative. They respectively either repel or attract a 
3-dimensional set of integral curves in the 5-dimensional phase space, the remaining two 
directions being degenerate. The points all correspond to f = singularities with the 
exception of the points with Yq = Vq which correspond to horizons. 

If g 2 > m + 1 or g 2 2 > m + 1 then in addition to the critical points for which Xq, Yq 
and Vo are all of the same sign, critical points of mixed signs also exist. These critical points 
are saddle points with respect to most trajectories in the phase space. The W = Z = 
solutions form one separatrix of trajectories with endpoints at the saddle points; for these 
solutions we still have r — > as the critical points are approached. Other separatrices are 
formed by the A = W = solutions, or the A = Z = solutions, as appropriate. For these 
solutions f — ■> oo as the saddle points are approached, as is discussed in Appendix B. For 
most trajectories, however, including the A > ones which are of prime interest to us, it is 
necessary to examine the behaviour of the solutions at the phase space infinity in order to 
determine the asymptotic (f — > oo) behaviour of solutions for which an asymptotic region 
exists. Apart from such solutions there are also many trajectories with two endpoints on 
the A = W = Z = curve which have no asymptotic region. Generally, they connect two 
points at which r — > 0. However, a subset with two endpoints in the V = Y , W = 'jZ 
subspace describe Schwarzschild-de Sitter like solutions, as was discussed in §2.2. 

All critical points at infinity other than the points L(y) and Si^ are found to corre- 
spond to f — > oo provided that < g 2 < m + 1 and < g 2 2 < m + 1. If g 1 = then 
the points P\,i correspond to r — > const, (indicating the presence of Robinson-Bertotti-like 
solutions), while the points AT 12 still correspond to r — > oo. The same is true for the points 
<5i,2 and Ri t 2 respectively if g 2 = 0. If > m + 1 then points AT 12 correspond to r — > 0, 
while if g 2 2 > m + 1 then the points similarly correspond to f — > 0. To discuss the 
asymptotic form of the solutions it is perhaps more convenient to examine the behaviour 
of the metric functions of the more usual Schwarzschild-type coordinates 

g ab dx a dx b = -e 2iL dt 2 + e 2i df 2 + ^g^dx^dx 13 , (2.43) 

rather than (2.1a). In Table 1 we display the asymptotic form of the metric functions 
(2.43), and of the scalar field, for integral curves from regions of the phase space at a finite 
distance from the origin which approach each of the critical points on the sphere at infinity 
for which f — > oo 

tt 8 If g 2 > m + 1 or g 2 2 > m + 1 then we will also have f — > oo for A = trajectories which approach 
points L(y) with y < 0. However, such trajectories are confined to the sphere at infinity and thus do not 
represent physical integral curves, and so we omit them. 
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Values of constants 


e 2n 


e 2 * 


e 2K<J 


M 1)2 


A > 


const. 


const. 


const. 


N lt2 


g x 2 < to + 1, A 1 > 0, A = 


f 2 


^(g^-1) 




Pit 


g 1 ^ 0, A 1 > 0, sign A = sign (g* - 1) 


f 2/ Sl 2 


const. 


fWgi 


Ql,2 


g 2 ^ 0, A 2 > 0, sign A = sign (g 2 2 - 1) 


f 2 / g 2 2 


const. 




R\,2 


g 2 2 < to + 1, A 2 > 0, A = 


f 2 


f 2(g 2 2 -l) 


f m S 2 


Tl,2 


A < 0, A = 


f 2 


f" 2 


const. 



Table 1 Asymptotic form of solutions for trajectories approaching critical points at 
infinity from within the sphere at infinity. 

In order to classify the various solutions we must first of all determine the nature of 
the various critical points at infinity. It is straightforward but laborious to evaluate the 
eigenvalue spectrum for small perturbations about the points. In Table 2 we summarise 
the results for such an analysis. We display the eigenvalues for the points with X > 0. For 
the corresponding points with X < the sign of the eigenvalues is simply reversed. 

The qualitative behaviour of the trajectories is largely dependent on the values of g 1 
and g 2 , apart from the case of points Mi p. These points - the only ones which correspond 
to solutions asymptotically flat in terms of f - are endpoints for a 3-dimensional set of 
solutions for all values of g 1 and g 2 . These solutions are just those lying in the W = Z = 
subspace, which is physically equivalent to A x = A 2 = . . . = A = 0. Thus models with 
non-zero A^ possess no solutions which are asymptotically flat in terms of f if A ^ 0. 

The eigenvalues for small perturbations near the points Nip, -Pi, 2, Qi,2 and Rip are 
essentially independent of the constants g v i > 3. The only exceptions is one eigenvalue 
at each point in the case that at least one of the a; or fa is less than one half. On account 
of the ordering (2.2) we also have 

< a 3 < a 4 < . . . < a s -i < a s < 1, (2.44) 

and 

1 > (3 3 > fa > . . . > Ps-i >(3 S >0. (2.45) 

Thus either ct 3 or (3 S has the smallest value of the and fa - this value being important 
in defining coordinates at the phase space infinity which lead to a well-defined spectrum 
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Eigenvalues (with degeneracies) 


m 


0, (2); 2; y; v.* 


Mi 


-i, (^); m 1 _ 1 , (*). 


Nt 


— (m + 1 -g 2 ) 


, (3); -Uf-l] 








Pi 


-i, (2); 4 


1 ± 


/9 + (m-9) gl 2 


. Si fei " 


-g ) 




Y 1 + (m - l) gl 2 


' l + (m- 


■l)g! 


2 ' 


Qi 


-i, (*); ^ 


1 ± 


/9 + (m - 9)g 2 2 " 


-g 2 (g/3 


-g 2< 


) 


Y 1 + (m - l)g 2 2 


1 + (m - 


■l)g 2 


2 ' 


Ri 


^( m+1 _ gj .), ( , )i i( gj »_ 1 ; 


) • ~ g2 ( 

' m V ^ 


-g 2 


)• 


Si 


-2; - 


1 J 


/ 8m( gl -g 2 )A 
-m± ]l m z — - 


1 




- (ro + 1) ,(*); = 


2 

> 









Table 2 Eigenvalues of critical points at infinity. The eigenvalues for small pertur- 
bations which are degenerate have the degeneracy listed in brackets. 
* The values of y and v listed are defined by (2.41) and V = vX in (2.42). 



of linearised perturbations. If we define 



, &3 , ifa 3 <ft, 
° Q 1 & otherwise, 



'2' 



and 



g s , if/3 s <« 3 , 

g 2 , otherwise, 1 ' ' 



then we find that one of the eigenvalues at the points AT 1]2 , Pi ,2 depends on the factor 
(gi ~~ g Q )) while one of the eigenvalues at the points <5i,2, FL\,2 depends on the factor 
(g^ — g 2 ). However, since both these factors are positive for each choice of a and (3 there 
is no qualitative difference between the alternatives. 
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In Table 3 we summarise the nature of the set of solutions with endpoints at iVi )2 
and -Pi, 2- The corresponding results for i? 1;2 and Qi )2 respectively may be obtained by 
substituting g 1 — > g 2 , — > Z. We display the dimension, <ix, of the maximal set, A, 
of trajectories with endpoints at each point - for g 2 < m + 1 this means the dimension 
of the set of trajectories attracted to (repelled from) the point Ni (AT 2 ), and vice- versa 
if g 2 > m + 1. In the case of the point Pi (P 2 ) it means the dimension of the set of 
trajectories attracted (repelled). Each point with dj, < 5 will also be the endpoint for a 
(5 — <i_4)-dimensional separatrix of saddle point trajectories. 



#1,2 


g x 2 <l 


Sl 2 = i 


1 < g x 2 < m + 1 


g x 2 > m + 1 




4 


4 


3 


5 


g 1 > Nature of 


W = 0, A > 0; 


W = 0, A > 0; 


W = 0, A = 


A > 0; 


solutions 




3-dim W=0, 






{A} 


W = 0, A < 


A < separatrix 




A < 




5 


5 


4 


4 


g 1 < Nature of 


A > 0; 


A > 0; 


A = 


W = 0, A > ; 


solutions 




4-dim A < 






{A} 


A < 


separatrix 




W = 0, A < 



Pl,2 


g x 2 <l 


g x 2 >l 


dA 


3 


4 


g 1 > Nature of 

solutions {^1} 


W = 0, A < 
separatrix 


W = 0, A>0 


dA 


4 


5 


g 1 < Nature of 

solutions {„4J 


A < 
separatrix 


A > 



Table 3 Nature of trajectories that approach points ATj.,2 and Pi,2- 

Of most interest are the points which are endpoints for a 5-dimensional set of trajectories, 
as they represent solutions with the most typical behaviour. If g 1 > then such points 
(AT lj2 ) exist only when g 2 > m + 1 and A 1 < 0. Such points exist for all g 1 < if A 1 > 0: 
for g i 2 < 1 they are the points iV 1)2 , and for g 2 > 1 the points Pi )2 . For the higher 
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derivative theories, however, g 1 > and g 2 < m + 1 and thus the points AT 1)2 or Pi 2 are 
endpoints for at most a 4-dimensional set of solutions. 

The nature of the set of solutions with endpoints at and is dependent on both 
the constants g. and the constants A.. However, the three eigenvalues corresponding to 
directions which lie within the (anti)-de Sitter subspace V = Y , W = 'jZ are independent 
of the g. and A., and so any differences are determined by the remaining two eigenvalues, 
which are given by the solutions of the equations 

A « +A * + 2(g 'mA S2)Ag = ° ^ 



for the points Si 2, and 



2 ,'m+l\ 2(m+l)(g 1 -g 2 )A g 



A T 2 + — A T + "/a " 9 = (2.486) 

for the points T 1)2 , where 

s s 

\ = ^ 2 g 2 7 2 + E \W 2ai = " E X iSAl 2ai ■ (2.48c) 

i=3 i=3 

For all choices we find only two possibilities: the dimension of the maximal set of trajec- 
tories with endpoints at Si ; 2 is either three or four; while the dimension of the maximal 
set of trajectories with endpoints at T12 is either four or five. The later was true in the 
case of the Kaluza-Klein models studied in II. 

This completes our classification of the solutions for Einstein gravity coupled to a 
scalar field with potential (1.10) for which the constant A, defined by (2.14), (2.15d) is 
non-zero. For non-zero X i solutions which are asymptotically flat in terms of f do not 
exist. This is not the case if A = 0, however, as we shall see in the next section. 



3. Solutions with a Schwarzschild subspace (Special polynomial R theories) 

If there exist solutions 7 of (2.14) such that the constant A defined by (2.15d) vanishes, 
then the structure of the phase space is significantly altered. No such solutions exist if 
s = 2. However, if s > 3, which applies in the case of the R + aR 2 theory and other 
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higher derivative models with potential (1.13), then such solutions are possible. We see 
from (2.17) that if A = we immediately obtain the Schwarzschild solution - thus all 
solutions lying in the 3-dimensional subspace V = Y , W = •yZ (for appropriate 7), and 
not just those in the Z = plane, are Schwarzschild solutions. This may be verified by 
direct integration since equations (2.15) are now equivalent to the equations which lead to 
the solutions of Appendix B.l, with the added restriction that c ± = mk. 

If we compare the phase space structure with that of the subspace of §2.2, we observe 
that the points <Si ; 2 and T12 do not exist in the Schwarzschild subspace. Instead, we have 
an additional one parameter family of critical points at a finite distance from the origin. 
These may be parameter ised in terms of their arbitrary value of Z = Zq > 0, and are 
located at 

X = Y = V = 0, Z = Z , W = <yZ . (3.1) 

We shall denote these points O(Z ), and the point at infinity (Z — > 00) will be denoted 
0\. These are the only additional critical points. 

Small perturbations about the points O(Z ) yield three zero eigenvalues, indicating 
a high degree of degeneracy. However, their properties with regard to the structure of 
the phase space may be ascertained from the following observations. Firstly, if Y = X 
or Y = (to + l)X/(m - 1), then P = 0, X' = and Y' = 0, so that the motion of the 
trajectories is entirely in the Z direction, as is depicted in Fig. 8. Furthermore, since (2.23) 
is now true for all solutions we may foliate the 3-dimensional subspace by a stack of planes 
Y = m(X + k)/{m — 1) to which trajectories are confined. Each plane is described by a 
2-dimensional autonomous system 

X' = X 2 -m 2 k 2 , (3.2a) 

771 

Z' = —— (X + k) Z, (3.26) 

TO — 1 

with solution 

Z = Co \X + mk\ 1/2 \X - mk\ {m+1)/[2(m - 1)] (3.3a) 

if k 7^ 0, or 

Z = C X m/ ^ m - l \ (3.36) 

if k = 0, where C is an arbitrary constant in both cases. The pattern of these trajectories 
is depicted in Fig. 9. The critical points at ±mk are of course those that lie on the 
A = W = Z = curve, and as expected they either attract or repel a 2-dimensional bunch 
of solutions in the subspace. From Figs. 8 and 9 we can see that for finite Z the points 
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O(Zq) neither attract nor repel any trajectories, and so there are no solutions for which 
they are endpoints. Clearly, however, all trajectories which reach infinity other than those 
in the planes depicted in Fig. 8 approach the point 0\. This is borne out by the plot of 
the hemisphere at infinity with coordinates 9 and <fi defined by (2.24), which is shown in 
Fig. 10(a). 

The fact that 0\ appears to be an endpoint for many solutions is in fact something 
of a misnomer, which arises from the fact that Z ~ X m 1 '( m_1 ) as X — > oo for all solutions 
(3.4) with an asymptotic region, so that Z grows more rapidly than either X or Y. Thus 
although the solutions (3.3) do not have X = 0or7 = 0asZ^ oo, they are nonetheless 
projected onto the north pole, Oi, if coordinates (2.24) are used. This degeneracy can be 
lifted if instead of (2.24) we use coordinates (p, 9, 4>) defined by 

X = psm9cos<i) Y = psm9sm<i) Z^ m ~ 1)/rn = pcos9 , (3.4) 

and perform the same analysis as before. This yields the plot Fig. 10(b) for the sphere at 
infinity. In these coordinates we obtain a line of critical points on the sphere at infinity 
with arbitrary angle 9, 

< 9 < |, 4> = arctan ^ m . (3.5) 

Perturbations about these critical points still yield one zero eigenvalue, corresponding to 
the degenerate direction 9, but each point in the first (third) quadrant is found to attract 
(repel) a 2-dimensional set of solutions from regions of the phase space at a finite distance 
from the origin. 

The crucial question now is: do any trajectories which lie outside the Schwarzschild 
subspace have an endpoint at 0\1 If such trajectories do exist, and if any of them curve 
back to another endpoint in the Schwarzschild subspace corresponding to an horizon, then 
we would have solutions with non-trivial scalar fields which violate the no-hair theorems. 
The corresponding Schwarzschild solution in any equivalent higher derivative theory would 
be non-unique. The answer is that the eigenvalues for perturbations in the two extra 
directions are given by 

A 2 = ^( gl -g 2 )A,, (3.6) 

where A g is given by (2.48c). Provided that A g > then 0\ will be a "centre" with 
respect to the extra directions, and the only solutions with an endpoint there will indeed 
be just the Schwarzschild solutions. If A < 0, however, then 0\ will be a saddle point 
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with respect to the extra two directions, and a 4-dimensional separatrix of solutions will 
have an endpoint at 0\ within the full 5-dimensional phase space. Unfortunately, it is 
not immediately obvious whether the solutions with an endpoint at 0\ which lie outside 
the Schwarzschild subspace have a second endpoint on the A = W = Z = curve at a 
point corresponding to an horizon or to a singularity. Thus we cannot extend the no-hair 
theorem, or the corresponding uniqueness theorem for higher derivative theories, to the 
case in which A g < 0. 

For the R + aR 2 theory, or indeed for any action (1.7) with coefficients (1.11), which 
gives rise to an equivalent potential (1.13), we find that a Schwarzschild subspace is ob- 
tained only for e = 1 and 7=1. Consequently, A g = m/ (8a\/m + 1) and 

2 — m 

X o = 4a(m + 1) ' (3 - 7) 

so that Oi is an endpoint for Schwarzschild solutions only provided a > 0. This is of 
course the same condition required for Whitt's proof in the case of the R + aR 2 theory 
(c.f. Appendix A). 

To complete our proof of the uniqueness theorem for higher derivative black holes it is 
still necessary to check the asymptotic behaviour of the metric functions near the critical 
points at infinity for which e 2Ka is not asymptotically constant in terms of the physical 
radial coordinate r instead of the coordinate f used in Table 1. The results for higher 
derivative actions with equivalent potential (1.13) are listed in Table 4, in terms of metric 
functions for the physical metric in Schwarzschild- type coordinates: 

exp ( ~ AKa ) g ah dx a dx b = -e 2u dt 2 + e 2v dr 2 + r 2 g afi dx a dx p . (3.8) 
\m^/m + 1/ 



(Near the points Ti 2 the solutions are of course still asymptotically anti-de Sitter in terms 
of the coordinate r as well as f.) It is evident that, as required, none of these solutions are 
asymptotically flat. 

From Table 3, we see that for the higher derivative theories points AT 1]2 are endpoints 
for the 3-dimensional set of W = 0, A = solutions, while points Pi ,2 are endpoints for 
the 4-dimensional set of W = 0, A > solutions. The constant g 2 can take either sign 
depending on the relative values of m and n, and its absolute value is only restricted by 
g 2 < m + 1. Thus the dimension of the maximal set of solutions with endpoints at Ri j2 
and Qi,2 is three, four or five depending on the particular values of m and n. For m = 2 
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Values of constants 


e 2u 


e 2v 


iv l,2 


prn — 1 

a 


r 2 




p 


pTL— 1 

< U, A > U 

a 


r 4/(m+2) 


const. 


<2l,2 


m^2(n-l), - <0, 
a 

sign A = sign [m — 4n(n — 1)] 


r 4(n-l)(2n-l)/[m-2(n-l)] 


const. 


Rl,2 


- < 0, A = 
a 


r 2 


r [m-4n(n-l)]/[(n-l)(2n-l)] 



Table 4 Asymptotic form of the higher derivative theory metric for polynomial R 
actions (1.7) with coefficients (1.11) in terms of the physical radial coordinate r, for 
the critical points at infinity at which r is not proportional to f. 

(i.e. the 4-dimensional theory) and n > 3, points R\^ are endpoints for a 5-dimensional set 
of solutions of either sign of A, while points <5i,2 are similarly endpoints for a 4-dimensional 
separatrix. If m = 2 and n = 2, the dimension of these sets is reduced by one with the 
additional restriction that Z = 0. In this case - or for any m and n such that m = 2(n — 1) 
- the points Ri$ correspond to asymptotically anti-de Sitter solutions. These solutions 
appear to form a special class unrelated to those of §2.2. 

The points Si^ and T 1>2 exist provided that m ^ 2(n — 1), 7 being given by roots of 
the polynomial (2.14): 

[2(n - 1) - m] (e7 2A ) " - 2n (e 7 2/fc ) " 1 + m + 2 = 0, (3.9) 

for which e^ 2 / n ^ 1, (and 7 > 0). If we factor out the Schwarzschild root (3.9) reduces to 
the polynomial 

. n-2 

[2(n - 1) - m] (e7 2/fc j - (m + 2) ^ (e~f 2/h ) = 0. (3.10) 

j=o 

For the R + aR 2 theory (n = 2), for example, the only solution (with 7 > 0) is 

* = -!, 7=g±§), ™>2, (3.11) 



- 32 - 



which gives 



7 - 1 A= 7 — ^ — , (3.12) 

' (m + 2){m - 2)a v y 



as is appropriate for the exact solution (2.19). Thus the points Si^ and Ti 2 exist only for 
a < 0, m > 2. We find that 

A = 7 m f ^1 (3.13) 

and consequently a 4-dimensional set of trajectories have endpoints at £1,2, and a 5- 
dimensional set of trajectories have endpoints at the points T12, with anti-de Sitter asymp- 
totics. If a > 0, the structure of the phase space is completely different: the points Si ,2 
and T12 do not exist. Instead we find a class of solutions with two endpoints lying in 
the V = Y, W = (m + 2)Z/(m — 2) subspace - these are asymptotically de-Sitter. In 
the case of four dimensions (m = 2) no subspaces of purely Schwarzschild-(anti) de Sitter 
solutions exist in the R + aR 2 theory. However, the points R\^ (which are endpoints for 
a 4-dimensional set of solutions) nonetheless have anti-de Sitter asymptotics. 



4. General polynomial R theories 

The analysis of the previous two sections applies exactly only to those polynomial 
R theories with coefficients (1.11). For lagrangians of degree n > 3 in R we shall see, 
however, that in the general case the field equations become equivalent to those given by 
one exponential sum potential or another near all the critical points. Consequently, the 
global properties of all theories polynomial in R can be obtained from the results of §2 and 
§3. As an independent check on this argument in one case, however, we shall first treat 
the general cubic theory in detail. 

4.1 The R + aR 2 + bR 3 theory 

The general cubic theory with f = R + aR 2 + bR 3 is described by a potential (1.6) 
which contains terms non-polynomial in exp [na/^/D — lj . The potential (1.6) may be 
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conveniently written as 



V = -r-^ < X 1 exp — + A 2 exp ' 



4k 2 



m 



m 



be A. 

M + AT exp 



3m 



+ A„ exp 



m 



where 



2ae ( 2 9 
A i = 276* 1° "2 6 



and £ is now defined by 



m + 2 

A, 



-(m - 4) 
— 2eea 



1 



\Jm + 1 



3^ 



a 



11/2 



A, 



1, if 1 + 3612/a > 0, 
-f, if 1 + 3612/a < 0. 



(4.1a) 
(4-16) 

36), 
(4.1c) 

(4.1d) 



As before, £ is defined by (1.2b). The quantities g 1? g 2 and g 3 are still identical to those 
given by (2.10a) for n = 3. However, the coefficients A^ now differ in general from those 
given by (2.10c) - the new coefficients A^ reduce to the values (2.10c), with A 4 = 0, in the 
special case a 2 = 36, when the potential is an exponential sum. Our new definition for e 
also reduces to the former definition (1.13b) if a 2 = 36. The Einstein-scalar field equations 
with the potential (4.1) are given by 



C" = (m - l) 2 Ae 2C + \e 211 + \e 



2X 



\be\ A 2 



3/2 



+ A 3 e 2 ^+ 2X )/ 3 , (4.2a) 



rj" = m (m - l)Ae 2C + — (m+l- g A A, e 2r] 

m V 1 / 



+ — (m+ 1 -g 1 g 2 ) A 2 e 
m 



-|(gi-g 2 )A 4 e 



2X 



\be\ A 2 



3/2 



i( m + i- gl g 3 )V^w 



2(2 V +X)/3 



\be\ A n 



1/2 



(4-26) 



x " = m{m - l)Ae 2C H (m + 1 - g 1 g 2 ) A^ 



+ ^(™ +1 -^)v 2X 



6e A 



4 p 2(r,-X) 



-i 3/2 



|6e| A 



+ ^(- + i-g 2 g 3 )V 2(r?+2X)/3 



^( gl -g 2 )A 4 e 2 ^/ 3 



N A 2 



1 1/2 



(4.2c) 



with the constraint 



(m + l)C ,2 + 



2mC(g 2V '-g lX ') ( 1 + (^~ 1 )g 2 2 ) /2 



Si §2 



-2 



(Si -g 2 ) 2 

(1 + (m - l) gl g 2 ) , , (l + (™ " ^Sx 2 ) 



(§i -S 2 ) 2 
A.e^ A.e 2X 



+ 



X ' 2 + (m-l)\e 



2C 



m m 



(gi " g 2 ) 2 



3/2 x 

_|_ ^3 e 2(r,+2X)/3 = Qj 

m 



(4.2d) 

where the variables £, ^ and x are defined as before by (2.3)-(2.5). If we now define the 
variables X, Y, Z : V and W by (2.7) once again, and use the constraint (4. 2d) to eliminate 
the e 2C - terms from (4.2a-c) we obtain the first-order system 



-1 



m 



x' = I (a^ 2 + A 9 r 2 + \w A ^z 2 /A - (m 1)p 

ml 1 z 6 ) m 

| ( gl 2 - l) A^ 2 + ( gl g 2 - 1) A 2 r 2 + ( gl g 3 - 1) A 3 ^/3 Z 2/3 

+ g 1 (g 1 -g 2 )A 4 ^ 4 / 3 r 2 / 3 |-p, 

V ' = ^{ (8 A - 1) \ Z2 + (g/ - V 2 + (g 2 g 3 - 1) A,U"' ! Z-- 

+ g 2 (g 1 -g 2 )A 4 z 4 / 3 r 2 / 3 l-p, 



(4.3a) 



(4.36) 



(4.3c) 



Z' = YZ, 

W = vw, 



(4.3d) 
(4.3e) 
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where 



r = 



J^VT^ + ^Z 4 / 3 
\be\ A n 



3/4 



(4-3/) 



and P is still given by (2.8f). 

The differences between equations (4.3) and (2.8) do not give rise to any significant 
changes to the analysis of §2 and §3. The position of all critical points W = and Z = 0, 
and in particular the W = 0, Z = X = surface, is unchanged, as are the eigenvalues for 
small perturbations about them. The number of critical points at infinity with Z ^ or 
W 7^ is also the same, and their location is more or less the same. For points and 
<5i )2 the definitions (2.39) and (2.40) are the same, if we now use A 2 as defined in (4.1c). 
These points now exist if ee/a < 0. In the case of points AT 1)2 and P\,i, we must make the 

1/2 

replacement X 1 — > X 1 + A 4 [A 4 /A 2 ] in the definitions (2.37) and (2.38), with X v A 2 and 
A 4 given by (4.2c). These points now exist if 



6 > 0, a 2 > 36, e = +1, 



ea < 0, 



(4.4) 



or if 



6>0, a 2 > 36, e=-l, (a 2 - 46) ea > 0. 



(4.5) 



In the case of points Si^, T\,i and Oi, our previous definitions are once again valid if we 
replace equations (2.14) and (2.15d), which respectively define 7 and A, by 



A lSl + A 2§2 



A 2 



-1 3/2 



+ A 3g 3 7 4/3 + A 4 ( gl -g 2 ) 



A 2 



nl/2 



= 0, 



(4.6) 



and 




i ±7 4/3 



3/2 



+ x^ 



(4.7) 



Here the plus (minus) sign corresponds to be > (be < 0). The eigenvalues about the 
various critical points are unchanged from those given in Table 2, if A is defined by 



7 



4/3 < A 



2 to 2 



1 



± 7 4 / 3 



1/2 2 1 

+ 3 A 3 g 3 + 3 A 4 (Si -g 2 ) 



A. 



± 7 4 / 3 



1-1/2- 



= - A iSi - 3 A 3 g 3 7 4/3 



± 7 4 /3 

A n 



-1/2 



^(2g 1 +g 2 )7 4/3 + ^g 1 



(4.8) 
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instead of (2.48c). 

A Schwarzschild subspace (A = 0) is obtained only for e = +1 and e = +1, the 
appropriate solution then being given by 7 = 1. Eigenvalues for small perturbations 
within the subspace give the same results as previously, and furthermore we find that for 
perturbations in the additional two directions the eigenvalues are once again 



— m 



(4.9) 



4a(m + 1) ' 

Thus just as in the special case of polynomial actions with coefficients (1.11), the black 
hole uniqueness theorem applies to solutions with a > 0. Most significantly, this result 
applies independently of the coefficient 6 of the R 3 term. 

With regard to the (anti)-de Sitter subspaces, we find that (4.6) has the solution 



(m + 2) \(m - 2)a 2 - 4(m - 4)6 ± a^/{m - 2) 2 a 2 - 4m(m - 4)6 



1 



4/3 



3(46 



2(m-4) 2 6e 



a 2 e 



if m 7^ 4, 

if m = 4, 
(4.10) 



where the sign of e is taken so as to make the r.h.s. positive. Consequently 



-(m - 2) a =F y/ (m - 2) 2 a 2 - 4m(m - 4)6 



7 " 4 / 3 A = 



-1 

3a 



2(m + 2)(m-4)6 



if m 7^ 4, 
if m = 4. 



(4.11) 



Thus exact Schwarzschild-de Sitter-type solutions of the form (2.18) are also obtained for 
m = 2, in contrast to the R+aR 2 theory. Furthermore, although in the case of the R+aR 2 
theory only anti-de Sitter solutions are admissable if a > 0, for the present theory both de 
Sitter and anti-de Sitter solutions are obtained for a > if m 7^ 4, (for either sign of 6). 



4.2 Theories of arbitrary degree 

At the next order with f = R + aR 2 + bR s + cR 4 the polynomial giving R in terms 
of /' has three branches, giving rise to potentials such as 



V = 



4k 2 c 3 



exp 



-2Dna 



(D-2)VD^1, 



256 



+ c 2 q - 26c 3 E Q 



(4.12a) 



6qc 



+ 3c 4 



+ 
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where 



q = 



64c 3 



16c 2 
Aabc - 6 3 



-ac 



8c eexp 



(. 



2ko 



and 



S ± = ( S o±\A 3 + S o 2 ) 



2/3 



(4.126) 
(4.12c) 

(4.12d) 



Such a potential would lead once again to equations similar to (4.3), if the g i of (2.10a) 
with n = 4 are used, but now with a further level of nesting of terms involving irrational 
roots. In the special case that 36 2 = 8ac, (i.e., q = 0), E + = 2E Q and £_ = 0, so that the 
potential (4.12) and the resulting differential equations become almost identical to those 
of the R + aR 2 + bR 3 theory, the only differences being the values of the A i and the factors 
a 3 and /3 3 , (which appear implicitly in (4.3) as the powers involved in the roots). 

We will not study the theory generated by the potential (4.12) in detail, but remark 
that it should not lead to any differences from our former analysis any more substantial 
than, say, the differences between the cubic and quadratic theories. In particular, the 
uniqueness of the Schwarzschild solutions should still apply to theories with a > 0. 

We will now show that the uniqueness theorem does indeed apply to theories with 
a polynomial R action of arbitrary degree n, if a 2 > 0. We note, first of all, that the 
Einstein-scalar field equations with a general potential V, derived from the action (1.4a) 
are given by 

" 2 (giX -g 2 7 7)' 



(" = ( m - l) 2 Ae 2C - 4 K 2 Vexp 
v " = m (m - l)Xe 2C - 



(m + l)AK 2 V dV 
m da 



x" = m(m — l)Ae 2 ^ 



(m+ 1)4k 2 V 



m 



dV 
da 



exp 


"2(g lX - 


-s 2 vY 


- Si- 


g 2 J 


exp 


"2( gl X- 






g 2 \ 



(4.13a) 
(4.136) 
(4.13c) 



with the constraint 



(m + l)C /2 + 



2mC'(g 2 ?/-g 1 X / ) , (l + (^~l)g 2 2 ) /2 o (1 + (m - l)g l8a ) , , 
+ ^ — 7 ^ — V -2 -, L r\X 



(gi " g 2 ) : 



(gi "g 2 ) ; 



(l + (m-l) gl 2 ) 



(gi -g 2 ) : 



X' 2 + (m-l)Ae 2C - 



4k 2 V 



exp 



2(g 1 X~g 2 ??) 

Si — s 2 



= 0, 
(4.13d) 
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where we have used the coordinates (2.1), and functions £, r\ and x defined by (2.3)-(2.5) 
with g 1 and g 2 as yet undetermined, (apart from the assumption that g 1 7^ g 2 ). In (4.13) 
a is assumed to be defined implicitly by the inverse relation to (2.4) and (2.5), viz. 

m f x-v \ u ^ 

Kff = y(iFi;J' (414) 

For a general polynomial R action (1.1), (1.7) this implicit definition becomes 
on account of (1.2), while 

v = h exp f 2 ^) 5> - 1)a ^ p ' (416) 

V 7 p=2 

and 

^ = ^ exp (z^)^ (gi _ pg3 ) bp, (4 . 17 ) 
da mn \ m I c " v 

if we take g v g 2 and g 3 to be given by (2.10a) as before. Thus if we define variables X, 
Y, Z, V and W by (2.7) we obtain the 5-dimensional autonomous system of first order 
differential equations 

X' = — Y(p - l)a R p - (m ~ 1)P , (4.18a) 

p=2 



k 



Y ' = ^r\ 81(81 - g 3 )^ + E [(s/ - 1) -p (s!g 3 - 1) 



4(m+ 1) 



p=2 

(m + 2)i2 + Yl t m + 2 ^ a p^ \ ~ P > 

p=2 

k 



a R p \ -P 

p 



(4.186) 



v ' = -^r\ 8 2 (si - s 3 )^ + E [(8182 - 1) - p (s 2 s 3 - 1)] %R p \ - p 

p=2 



-eZ 



4(fc _l )(m + l) I ( 2 ( fc - 1) - m) R + E [~( 2k + ™) + 2 ( 2fc " j " P > 

(4.18c) 
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Z' = YZ, 

w = vw, 

where R is now denned implicitly in terms of W and Z by 

2(n-l)/n 

= 1 + 

p=2 



(4.18d) 
(4.18e) 



W 
~Z 



(4.18/) 



and P is still denned by (2.8f). 

The properties of the phase space defined by (4.18) differ little in their essentials from 
the examples studied earlier. We first note that an (anti)-de Sitter subspace is defined by 
the (constant) values of R given by roots of the polynomial 



k 

m + ^2[m-2(p- 1)] a pJ R p_1 = 0. 

p=2 



(4.19) 



The 3-dimensional subspace is given by V = Y and W = 7Z, where by (4.18f) 



e 7 2(n - 1)/n = l + E^ P ' 

p=2 



(4.20) 



R being a solution of (4.19). The field equations of the subspace are once again given by 
(2.15) with A now defined by 



A = -}Z(p- 1 K RP - 

p=2 



(4.21) 



Our earlier results concerning the (anti)-de Sitter subspace follow through. In particular, 
we are lead to the Schwarzschild-de Sitter type solutions (2.18), which in terms of the 
original physical metric are given by 



ds 2 = -Adt 2 + A _1 dr 2 + r 2 ^dx Q d^, 



with 



A = A 



2GM 


Ep= 2 b-i)vR p 




r m — l , , 

' m[m + 1) 





(4.22a) 



(4.226) 
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R being a solution of (4.19). Furthermore, we once again obtain Robinson-Bertotti solu- 
tions of the form (2.20) for our new definitions of 7 and A. 

Let us now consider the critical points of the system (4.18). It can be readily seen 
from (4.18a-c) that the only critical points at finite values of X, Y, V, Z and W must 
all have: P = 0, and also either (i) W = Z = H9 ; or (ii) R + J2p^P a p RP = and 
R + Xlp=2 a p R P = 0- The oru y value of R which simultaneously satisfies both of the 
conditions (ii) is R = 0. We therefore retrieve the same critical points as were found in §2 
and §3, viz. 

(i) . The familiar W = 0, Z = 0, A = surface discussed in §2. 

(ii) . X = Y = V = 0, W = Z = Z , (R = 0). These are the points O(Z ) which lie in the 
Schwarzschild subspace, discussed in §3. 

The phase space at infinity may be studied either by a direct analysis of the 4-sphere 
at infinity, using an approach similar to that described in §2.2 and §2.3, or by the compu- 
tationally more simple method of defining new variables 5, y, v, z and w by 

and classifying the 6 = critical points of the resulting field equations 

dr m L J 



±— = - [(y - l): 2 ./, - ,V-'-'./ 2 + «(("» - 1)1/ - m )] , (4-246) 



dy 

dr m 



±? = - [(« " i)^/! " g 2 ^/ 2 + «((rr» - 1)« - m)] , (4.24c) 
dr m 

±p- = - [my + z 2 f, + (m - l)g] , (4.24rf) 
dr m L 1 J 

±— = — \mv + z 2 f, +(m-l)q] , (4.24e) 
dr m L 1 J 



^ 9 Ostensibly it would seem that we can simply take Z = here. However, since ii is defined implicity in 
terms of W/Z the summation terms in (4.18a-c) will not vanish simply if the overall factor of Z 2 vanishes. 
The leading order R n term within the summations is of order W 2 Z~ 2 , and so for consistency one must 
also require that W = when setting Z = 0. 
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where dr = 5 1 d£, g = 5 2 P and 



A^e^Cp-lJa^, (4.24/) 

p=2 



/ 2 = £ 



(g 1 -g3) J R + ^(g 1 -pg 3 ) aW 

P=2 



{A.2Ag) 



with i? now defined implicity in terms of w and z. This method picks out all critical points 
at infinity apart from those with X = 0. To obtain all critical points we must repeat the 
calculation with each of the phase space coordinates in turn defined as ±1/5. 

It is convenient to divide the critical points at infinity arising from this analysis into 
three classes according to the behaviour of the scalar curvature R at the points: (i) R = 0; 
(ii) R = const. ^ 0; or (iii) R — > oo. 

(i) . Critical points with R = have f 1 = and f 2 = 0. Consequently the only possibilities 
admitted by the equations (4.24) are: the points L(y), defined by (2.41), (2.42); and the 
points Mi ; 2, defined by (2.29) (with V = Y and W = also). If we repeat the above 
analysis for the system of equations defined by putting Z = 1/5, with the other variables 
proportional to Z, then for R — we also obtain the point Oi defined by the Z — > oo 
limit of (3.1). 

(ii) . If R = const. ^ then we must have z ^ and to ^ on account of (4.18f). 
Equations (4.24) then imply that v = y and f 2 = 0. This latter condition is equivalent to 
(4.19) for R 7^ 0. Thus such critical points must lie in the anti-de Sitter subspace, which 
as we have already discussed is well-defined. By the analysis of §2.2, these critical points 
are therefore Si^ and Ti j2 , where the 7 and A of the defining relation (2.30b) and (2.31b) 
are now taken to be given by (4.20) and (4.21) respectively. 

(iii) . Obviously points at which R — > 00 cannot correspond to regular solutions. However, 
it is nevertheless interesting to check whether the structure of the phase space is preserved 
when compared to our earlier examples. If R — > 00, then on account of (4.15) 



R~± 



1 



na 

n 



e exp - 1 



l/(n-l) 

(4.25) 



In the limit R — > 00 the field equations are therefore equivalent to those derived from an 
exponential sum potential (1.10) consisting of a n + 1 terms, with g v g 2 and g 3 given by 
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(2.10a) and 

2(n - 1) + m(j - 3) A tl . 

g 3 -= L ^ /= =r> 4<j<n+l. (4.26) 
J 2(n — l)\/m + 1 

We could also derive the values of the appropriate coefficients, A., however, they are not 
of much concern to us here. We merely note that since the field equations are equivalent 
to those of an exponential sum potential, the only critical points at which R — > oo are 
precisely the points N\^, Pi, 2, Qi,2 and Ri,2, (with appropriate A. in the definitions). The 
asymptotic form of the solutions is therefore given by Table 1, or in terms of the physical 
metric, by Table 4. In particular, all of these critical points correspond to asymptotic 
regions in which the physical metric is not asymptotically flat. 

Thus the only integral curves which join a critical point corresponding to an horizon 
to a critical point corresponding to an asymptotically flat region must be the trajectories 
with one endpoint on the curve formed by the intersection of the W = 0, Z = 0, A = 
surface with the Schwarzschild subspace, and with a second endpoint at or 0\. It 

therefore only remains to determine the dimension of the set of such solutions. As before, 
this may be found by a linearised analysis of small perturbations about the points. This 
analysis will be unchanged from that of §2 and §3 since by (4.15) we have to leading order 
in R 

2k<j 



\Jm + 1 



l + 2aR + 0(R z ). (4.27) 



Consequently, the linearised perturbation equations will be identical to those obtain in the 
R + aR 2 theory. In particular, if a 2 > then no trajectories in directions orthogonal to the 
Schwarzschild subspace will have endpoints at or 0\. This completes our proof of the 
uniqueness theorem. We have of course assumed that a 2 7^ throughout. If a 2 = then 
it seems likely that the properties of the asymptotically flat solutions will be essentially 
determined by the value of a , where p is the least value such that a p 7^ 0. 



5. Discussion 

To conclude, we have shown that in theories with an action polynomial in the Ricci 
scalar, the only static spherically symmetric solution with a regular horizon is the Schwarz- 
schild solution, provided that the coefficient a 2 of the quadratic term is positive. In fact, 
if we drop the condition of regularity on the horizon then the only asymptotically flat 



solutions are still the positive and negative mass Schwarzschild solutions. The generalised 
scalar potential model of §2 does possess non- Schwarzschild W = 0, Z = 0, A > solutions 
with naked singularities which approach flat space asymptotically near the points M\$. 
However, these solutions correspond physically to A 1 = A 2 = . . . = A s = 0, a possibil- 
ity which is not admitted in the higher derivative theories. We have not considered the 
theories for which the quadratic terms vanishes but which have non-zero terms at higher 
order. In such cases the condition a 2 > presumably translates into a condition on the 
coefficients of higher terms in the series. The question of solutions with a 2 < is also not 
fully clear, and could perhaps be resolved by a numerical study. 

The fact that we have been able to prove a uniqueness theorem for static spherically 
symmetric black holes in the case that a 2 > has immediate important physical impli- 
cations when considered in conjunction with the earlier work of Pechlaner and Sexl [19] 
and Michel [20] . These authors observed that in the case of the R + aR 2 theory in four 
dimensions the solutions with R = are not the solutions which correspond to the weak 
field limit about any physical body such as a star or point particle. This is also the case for 
the general theory with a polynomial R action which we are considering here. Suppose, for 
example, that we wish to match our solutions onto a star with energy-momentum tensor 
T a b in its interior. If we add such a term to (1.8) and trace the result we find that 

^(D - 2)R + % { (p ~ f ) R p + (D- l)p(p - l)R p ~ 3 (RDR+ (p - 2)R C R, C 

= 2k 2 T. 

(5.1) 

Thus if we set T = at the surface of the star we cannot conclude, as in the case of the 
Schwarzschild solution, that R = 0. Our results show, however, that if a 2 > and R ^ 
anywhere in the domain of outer communications, then the solutions are not asymptotically 
flat. 

Pechlaner and Sexl, on the contrary, assumed the existence of non-Schwarzschild so- 
lutions which are asymptotically flat. Our full non-linear analysis seems to invalidate this 
assumption in the a 2 > case. Consequently, their weak field analysis of the R + aR 2 
theory in four dimensions, in which they derived experimental bounds on the parameter a 
from resulting fifth-force type effects, needs to be reexamined with appropriately changed 
boundary conditions. 

If the Schwarzschild solutions are not the solutions of physical interest in these models, 
it would seem that the large class of asymptotically anti-de Sitter and de Sitter solutions 
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hold more promise. These include both the exact solutions (2.18), (with A given, for 
example, by (3.12) and (4.11) in the quadratic and cubic order theories respectively), 
and also other solutions asymptotic to them at infinity, or at the de Sitter cosmological 
horizon. Asymptotically (anti)-de Sitter solutions have also been found in a number of 
models in D > 4 dimensions which incorporate a Gauss-Bonnet term [32-34] and other 
dimensionally continued Euler densities [35], and thus appear to be a generic feature of 
higher derivative theories. In fact, maximally symmetric solutions have been found to 
exist in a much wider class of higher derivative models [31], so presumably such models 
also possess asymptotically (anti)-de Sitter black hole solutions. 

Returning to the models with dimensionally continued Euler densities, we note that 
similarly to the solutions found here, only asymptotically anti-de Sitter solutions are found 
if the coefficient of the Gauss-Bonnet term is positive in the quadratic theory, while a de 
Sitter branch can be obtained at higher order [35]. One important difference between our 
solutions and those of Boulware and Deser [32], for example, is that the asymptotically 
anti-de Sitter branch of the Boulware-Deser solutions has a negative gravitational mass, 
giving rise to an instability. Our solutions have a positive gravitational mass, and so the 
question of their stability is still an open problem. The issue of the weak field limit has 
not been addressed in refs. [32-34], since of course it is really only relevant in compactified 
models. 

If the de Sitter sector is to be treated as a serious physical model it is clear that the 
effective cosmological term must be small to be consistent with observation. It is interesting 
to note that recent astronomical evidence actually favours a small positive cosmological 
constant^ 10 [36]. If higher order terms in R are obtained from the dimensional reduction of 
an action corresponding to the low energy limit of a higher-dimensional string theory, then 
the status of the effective cosmological term is at best uncertain. If we ignore the values of 
the dilaton and compacton and assume that the compact ificat ion scale is approximately 
of order \fa i (in units in which c = % = 1, as used throughout this paper), where a' 
is the Regge slope parameter, then the coefficients a p are of order (^pi anc k) 2 ^ p-1 \ an d 
A is of order (^ Planck ) -2 , to within a few orders of magnitude. Such a colossally large 
effective cosmological term would of course spell disaster for these models. However, no 
definite statements can be made without some knowledge of the expectation values of 
the dilaton and compacton fields, both of which couple non-trivially to the higher order 

* 10 The "best-fit" value favoured in ref. [36] is A = 3.1 x 1CT 52 to- 2 , or A = 8.1 x 1CT 122 in dimensionless 
Planck units. 
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curvature terms in four dimensions. In fact, these scalar fields should really be treated 
dynamically, which would necessitate a complete reexamination of the model. In the D > 4 
uncompactified quadratic order theory with a Gauss-Bonnet term the dilaton has the effect 
of removing the anti-de Sitter branch [34] , but it is not clear how the dilaton would affect 
compactified models. 

In view of these problems it would also be interesting, especially in the quadratic order 
case in four dimensions, to consider the addition to the lagrangian of a term comprising 
the square of the Ricci tensor. In that case the effective theory contains a massive spin 
2 field with a non-trivial coupling to gravity in addition to the scalar field. The effective 
energy momentum tensor of the extra excitations does not satisfy criteria usually required 
in the proof of the no-hair theorems, and thus it seems plausible that the no-hair theorems 
could be circumvented in such a model. However, the dynamical system arising from a 
static spherically symmetric ansatz for the metric and the other fields is considerably more 
complicated than in the models we have studied in this paper, and it is not clear to us 
whether it can be reduced to a tractable form. 

We remark parenthetically that the results of this paper seem to hold up some hope for 
the problem of finding black hole solutions in dimensionally reduced theories in which the 
internal space is non-Ricci flat, since such models can also be treated by the formalism of §2 
and §3. In I and II the potential corresponding to (1.10) was limited to two terms at most. 
However, we have seen that at least three exponential terms are required in the potential 
in order to obtain a Schwarzschild subspace, and thus no asymptotically flat solutions were 
found in I and II for internal spaces of non-zero curvature. More complicated models with 
s > 3, (or some equivalent condition if more than one scalar field is present), may yield 
more interesting results. 



Appendix A 

In the case of the quadratic action R + aR 2 in four dimensions with a > 0, Whitt [21] 
has shown that provided the energy momentum of any additional matter fields satisfies 
certain conditions then the usual no hair theorem for stationary, axisymmetric, asymptot- 
ically flat black holes [37] is still valid. He was able to obtain this result by proving that 
the curvature scalar for such solutions vanishes in the domain of outer communications, 
and therefore the effective theory coincides with the usual Einstein-Hilbert theory. Unfor- 
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tunately his proof cannot be generalised to more general polynomial actions of the form 
(1.7). For completeness we will present an account of his proof here, in order to show why 
the theorem cannot be directly extended to more general polynomial R actions, and also 
to correct some errors present in the original paper. 

Whitt's proof is based on a study of the 4-dimensional Killing bivector 



p ab = 2rri[ a k b ], (A.l) 

where m a is the spacelike Killing vector with period 2ir and k a is the timelike Killing vector 
normalized to unity at infinity. These vector fields commute [39] : 

k b m a . b - m b k a . b = 0. (A2) 

We will take l a to denote either of the Killing vectors k a and m a . Making use of (A. 2) and 
the identity 

2 

(^[a;^c]) ;C = ^ C Rc[Jb] , (A3) 

which is valid for any Killing vector [40], we find 

(ha;bk c m d] y d = -h d R d[a k b m c] . (A4) 

We can now use the vacuum field equations of the higher derivative theory to evaluate the 
right hand side of this expression. For the R + aR 2 case one finds 

(h^ m ^= - al [^ m \ (AM, 
or in the general case, with field equations (1.8), 

. -(ELv^;*^]) 

(l[a;bk C m d] y d = — V T V 1 ' ( A - 5b ) 

2 (l + Ep =2 a pP RP-i) 

Using the orthogonality of the Killing vectors to R- a and the antisymmetry of the deriva- 
tives of the Killing vectors, these expressions become respectively 

{l[a;bk c m d] y = 1 + 2a R ^ ' 
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and 

— Y] m 9 a p(R p ~ 1 )' ,d h t ]. n kb'm r ] 
(ha; b k c m d] y d = ^ P ~ 2 pP \ C] (AM) 

2(i + EU a P P RP ~ 1 ) 

The solution of equations (A. 6) is given respectively by 

l[a;bk c m d] = C(l + 2aR)- 1/2 e abcd (A.7a) 

and 

/ k N-V2 

l[a-bk c m d] = C I 1 + J^cyLR p 1 J e a&cd (A76) 

where C is an arbitrary constant. Now we know that the left hand side of (A. 7) vanishes 
on the axis of rotation since m a = there. On account of the asymptotic flatness of the 
solutions R cannot be singular everywhere on the axis, and hence C must vanish giving 

k[ a -,bPcd] = 0, m [a . b p cd] = 0, (A8) 

throughout the domain of outer communications. 

From (A. 8) it follows [40] that if the domain of outer communications is simply con- 
nected and admits no closed timelike curves, then p ab is timelike throughout the domain of 
outer communications, becoming null on its inner boundary, which is a null hypersurface. 

At this point Whitt [21] considers the trace of the field equations in the absence of 
matter, which for the quadratic theory yield 

-6a.R ;a . a + R = 0. ( A9) 

One can now multiply by R and integrate over the domain of outer communications to 
obtain 

6a J RR. a dE a = 6a J R> a R ]a y/=gd*x + J R 2 v ^gd 4 x. (AlO) 

The left hand side vanishes because R- a is orthogonal to the inner boundary and is zero on 
the outer boundary due to asymptotic flatness. Moreover the right hand side is positive 
definite, provided a > 0, since R. a cannot be timelike anywhere in the domain of outer 
communications by virtue of being orthogonal to p a b- Hence R must be identically zero 
there. 
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The conformally rescaled theory coincides therefore with the usual Einstein theory 
since the scalar field a is clearly zero if R vanishes and it follows that the uniqueness the- 
orem for stationary, axisymmetric, asymptotically flat solutions still holds. Moreover, it is 
easy to see that the result remains valid when additional matter fields are included in the 
action, provided that the energy momentum tensor is traceless, T a a = 0, and satisfies the 
matter circularity condition [40] T d [ a fc&m c ] = 0. In such cases the usual no-hair theorems 
are still valid. These conditions are always satisfied, in particular, by a stationary, axisym- 
metric electromagnetic field. Unfortunately, however, Whitt's result cannot be extended 
to the general polynomial R action since in that case the argument for R being zero in the 
domain of outer communications breaks down. Specifically, the trace of the field equation 
now becomes 

k 

-R + J2%[(p-2)R p + MR p - 1 y a , a ]=o (Ail) 

and therefore, integrating by parts and discarding the boundary term as before we find 
# 2 v/^d 4 x + 5>-2)a p / ^ +1 v^d 4 x = ^3p(p-l)a p / R^ 2 R, a R a ^d'x 

p=2 J p=2 J 

(A12) 

The integrals in the sum on the left hand side do not have a definite sign for even p, while 
the integrals in the sum on the right hand side do not have a definite sign for odd p. Thus 
it is impossible to conclude from this relation that R must vanish in the general case, or 
indeed in any special case other than n = 2, a 2 > 0. 



Appendix B 

We list here the exact solutions obtained in the cases (i) V = 0, i.e., A x = A 2 = 
. . . A = - appropriate to the W = 0, Z = subspace; (ii) s = 1, A = - appropriate to 
the A = surface in the W = subspace. Corresponding solutions for the A = surface 
in the Z = subspace may be obtained by substituting g-,^ — > g 2 , A x — > A 2 . 

B.l Solutions with W = and Z = 

The one-parameter family of solutions was derived in I. In that case solutions were 
parametrised in terms of the number n e of extra dimensions in the Kaluza-Klein model. 
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Here a parametrisation in terms of g 1 (or g 2 ) is more natural. We will use the former 
parametrisation. We find 



CAi 

^7 



exp 



1 



(m - 1)C£ 



m 



-(C + H) + const. 



[m 



m — 1 

- l)CAi exp [i(m - 1)C£] 



(m - l)^ 2 exp [(m - 1)CC] + (a Q + \{m - 1)C) Ai 



r = — — : — exp 



A A 3 



1 



a + -(m-l)CH 



2g 1 na (m - l)g 1 2 (c 1 - mfc) 



1 + (m - l)g 1 2 



+ const., 



where 



Ai = A - A x exp [(m - 1)C(] , 
while fc, A Q and Ai are arbitrary constants, C is a non-zero constant given by 

- A {m - If (l + (m - l) gl 2 ) C 2 = mV + (m - l)^ 2 , 
and the constant a Q is defined by 

mk + (m — l)g l 2 c 1 
tt ° = 1 + (m - l) gl 2 ' 



(5.1) 

(5.2) 
(5.3) 
(5.4) 

(5.5) 

(5.6) 

(5.7) 
(5.8) 



(5.9) 



and lies in the range — |(m — 1)|C| < a Q < \{m — 1)|C|. Solutions have an asymptotic 
(f — > oo) region only if A > 0. 

As £ — > — oo we find f — > 0, giving rise to a singularity except in the special case 
c 1 = mk = a Q = — \{m — 1)|C|, when f — > const., suggesting the possible presence of an 
horizon. This can be immediately verified to be true since (B.5) can be inverted. If we 
choose* 111 

M^ 1 ' c<°o: (Bio > 



" 11 Equations (B.10)-(B.12) correct a sign discrepancy in I for the A < and £ — > +oo cases. 
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we find 

e " = e " 25 = 1 - p|M • (B ' n) 

which correspond to the domain of outer communications of the positive mass Schwarz- 
schild solution for A > 0. Similarly, as £ — > +oo we find f — > 0, giving rise to a singularity, 
unless = m/c = a Q = |(m — 1)|C|. In that case we once again obtain (B.ll) if we choose 



B.2 Solutions with W = and A = 



These solutions may be derived as in I. If g 2 < m + 1 we find 



exp 



(5.13) 



exp 



777.(77 + ^0 
777 + 1 - g x 2 



-C + mk- (m- 1)6 £ 



777 



e = 



+ 1 - gl 2 ) So {a^ (CSi exp [(|(7 + fc) £])' 



l/(m+l- gl 2 ) 



l/(m+l- gl 2 ) 



1 \CB 1 
t = — — { — — exp 



B0B2 1 A2 



'2 gl na\ 1 CB 1 

exp \ = — < — — exp 

ml B 2 A 2 



2^ + 6 oK 



l/(m+l- gl 2 ) 



gl 2 /(m+l- gl 2 ) 



where 



with 



A ^TO> +1 - g i 2 K 



(5.14) 
(5.15) 

(5.16) 
(5.17) 
(5.18) 

(5.19a) 
(5.196) 
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while B Q , Bi, B 2 and k are arbitrary constants, C is a non-zero constant given by 

+ l) gl 2 ) C 2 = m 2 k 2 + (m + 1 - g 2 ) g 2 c 2 , (B.20) 
and the constant b is defined by 

(g 2 - l)mk + (m + 1 - g 2 ) g 2 c 
b n = — -A , = 7 . (S.21) 

and lies in the range — \\C\ < b Q < \\C\. 

We now find that the limit £ — > — oo corresponds to f — > except in the special 
instances when mk/(g 2 — 1) = c 1 = 6 Q = — 1|(7|, for which f — * const., suggesting the 
possible presence of an horizon. This indeed the case: (B.17) can be inverted and if we 
make the choice 

D. J -"0 2 ) O > U, /n 9 ",\ 

jyi "\A 1 /( J B m+1 S 2 m ), (7<0, ' 

then we find the solution 

ds 2 = -f 2 Mt 2 + f 2(g i 2 - 1} i- + f 2 g a0 dx a dx p , (S.24a) 



where 



A = A I 1 ^3 I ( S .246) 

lA^m + l-g^f^ 1 ^! 2 ' 

while the scalar field is given by 

2 

exp = ? . S-24c 

\ m / m + 1 — g z 

Similarly the limit £ — > +oo also corresponds to r — > except for the special cases when 
mk/(g 2 — 1) = c 1 = 6 Q = which f — > const. Equation (B.17) can once again be 

inverted, and we retrieve the solution (2.18) if we now make the choice (B.22) and* 112 

_ f -kJ{B™^B C > 0, 
Sl " 1 -B^Bf, C < 0, ' (S - 25) 



J12 There are some sign discrepancies in I and II for the expressions corresponding to (B.23), (B.24b) 
and (B.25). 
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The spacetimes thus have naked singularities except in the special cases above. 

An asymptotic region is defined only for A^^ > 0: f — > oo when (7£ = In \A 1 /B±\. This 
limit is approached at the points iV~i 5 2 at infinity (c.f. (2.37)). All A 1 > solutions have 
one endpoint at N± or N2, and another endpoint on the A 1 = curve of critical points. 
The A 1 < solutions have one endpoint on the A 1 = curve in the first quadrant, and 
another endpoint on the A 1 = curve in the third quadrant. The asymptotic form of all 
the X 1 > solutions (B.15)-(B.18) is given by 

e 2ti ~f 2 , e M ~f 2(g i 2_1) , e 2K<J ~f mg i. (S.26) 

Thus none of the solutions is asymptotically flat. The general solution holds for C ^ 0. 
If one sets (7 = while integrating the differential equations one finds a solution which 
corresponds to the (7 = limit of (B.24). 

If g 2 > m + 1 then the solutions given by expressions (B.13)-(B.21) are still valid 
but b Q now lies in the range b Q < — 1](7|, b Q > ||(7|. Furthermore, the behaviour of the 
solutions and the direction of trajectories near the critical points is greatly changed in 
some instances. We now find that the limit £ — > — 00 corresponds to f — > as before if 
b Q < — \\C\. However, if b Q > \\C\ we find that f — > 00. Similarly, the limit £ — > +00 
corresponds to f — > if b Q > \\C\, and r — >ooif6 Q <— \\C\. Critical points in the first 
and third quadrants have r — > 0, while those in the fourth and second quadrants have 
f — > 00. Since each solution has a different endpoint in the f — > 00 regime, the asymptotic 
behaviours vary. We find 

where 



±11 


C| + mk - 


(m-l)b Q 


^o±i 


C\ 





% = 21 jt-tt^, ^, (5.276) 



m(/c±||C|) 

*^WT' (B ' 27C) 

±l\C\ + ±{mk + b ) 

°° = 6~±~I|<3j ■ (B ' 27d) 
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and the upper (lower) sign refers to £ — > — oo (£ — > +00). For the limiting case solutions 
with 6 Q = ^\C\, if £ — > —00, or b Q = — 1|(7|, if £ — > +00 we obtain the asymptotic behaviour 
appropriate to (B.24), viz. 



(5.28) 



The limit, C£ = ln\A 1 /B 1 2 \ for the A x < solutions, which is reached at the points Ni^, 
now corresponds to f — > 0. The A 1 < solutions given by (B.13)-(B.21) all have one 
endpoint at Ni or N2 and another endpoint on the A 1 = curve. If this second endpoint is 
in the fourth or second quadrant the solutions have an asymptotic region, otherwise they 
do not. The A 1 > solutions, on the other hand, all have asymptotic regions: they now 
describe trajectories with one endpoint on the A = curve in the first (or third) quadrant, 
and a second endpoint on the same curve in the fourth (or second) quadrant. 

In addition to the above solutions, a second group of solutions also exist if g 2 > m+1 
and A 1 < 0. These solutions are given by 



\C\ 



= B exp 



-m(r) + fcg 
g 2 — m — 1 



(5.29) 



(5.30) 



= 5 m 5 2 (m - 



" 1} {A^ICT 1 cos - Co)) exp [((to - 1)6 - mk) £] } gl 



l/( gl 2 -m-l) 



(5.31) 



e v = 



-(g 2 -m-l)5 , . v _ . m /( g 2 - m -i) 

^A 1 1 / 2 |C|- 1 cos(c(^-^ ))exp(-^)} 1 



6 +(7tan ((7(£ - £ ) 

f = 5^ K 172 !^ -1 cos " ex P(- & oO} 1/(Sl2 " m_1) 

ex p(^^) = ^{[ A 1 1/2 ic , r 1 cos((5(e-g)] Sl2 exp[-(mfc+6 )e]} MBl 



(5.32) 
(5.33) 

l/(g 2 -m-l) 



(5.34) 

where A 1 and b Q are defined as before, £ , 5o, 52 and k are arbitrary constants, and C is 
a non-zero constant given by 



(l + (m - l)g x 2 ) C 2 = (g 2 - to - l) gl 2 Cl 2 - m 2 k 2 . 



(5.35) 
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These solutions have no asymptotic region, and correspond to trajectories with one end- 
point at Ni and a second endpoint at N 2 . These endpoints are reached when £ = 
£ + (2n + l)/(2ir\C\), for integer n. 



If g 2 = m + 1 we find the solutions 

e" = B 1 C 2 e^ 2 , 
e c = S exp k£ + X^e^ 

= 5 m 5 2 (m_1) jexp [(jfe - (m - 1)& ) £ + A 1 5 i e^j } 



l/m 



mBo exp 



k + b Q + X 1 B 1 Ce c ^ 



B0B2 



exp 



exp 



2g 1 ^o- 
m 



Bo 



exp 



l/m 



(m+l)/m 



(5.36) 
(5.37) 

(5.38) 
(5.39) 
(5.40) 

(5.41) 



where 5 , 5i 5 2 are arbitrary constants, 



>o = M§ + (™ + iK|> 



(5.42) 



fe is a constant which lies in the range |fc| > V m 2 — 1 joj/m, and C is a non-zero constant 
given by 

2 



C 



m — 1 



mk ± Jm 2 k 2 — (m 2 — l)c 2 



(5.43) 



and has the same sign as k. We now find that for /c > as £ ^ — oof^O except in 
the case that k = —c 1 = —b Q = (7/2, when f — > const. The same is true for fe < in the 
limit £ — > +00. On the other hand, if k < and £ — > —00, or if > and £ — > +00, 
then f — > 00. All solutions have an asymptotic region. In terms of the phase space, all 
trajectories approach the points iV lj2 at infinity which coincide with points L 5 j. In the 
case that k = —c x = —b Q = C/2 we can invert (B.40), and if we make the choice 



5 9 = 5 



> 



B - mB ° 2 
1 X ± C 



(5.44) 
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we obtain the solution 



ds 2 = -r 2 C In rdt 2 + 



ty> 2 771 J ,y, 2 



+ f 2 g a pdx a dx p , 



(5.45a) 



Clnf 



with scalar field 



exp 




,m+l 



(5.456) 
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Figure captions 

Fig. 1 Trajectories in the W = 0, Z = 0, V = Y plane. The bold lines Y = X and 

Y = (m + l)X/(m — 1) represent critical points which respectively correspond to the 
horizons of the positive mass Schwarzschild solutions (dashed lines), and the singularities 
of the negative mass Schwarzschild solutions. The trajectory through the origin is flat 
space. 

Fig. 2 The plane V = Y = X, W = 7Z, with radial coordinate p and angular coordinate 
8: (a) A < 0, A < 0; (b) A > 0, A > 0. The solutions are given in (2.20) in terms of 
Z = pcos 8/(1 - p). 

Fig. 3 The hemisphere at infinity for the 3-dimensional (anti)-de Sitter subspace: (a) 
A < 0; (b) A > 0. 

Fig. 4 The projection of trajectories in the W = 0, Z = subspace onto the X,Y plane, 
with V given by (2.12), for non-zero c x : (a) g 2 < m + 1; (b) g l 2 > m + 1. The broken 
line corresponds to the Schwarzschild solution. The bold lines represent sections through 
the cone of critical points A = 0, W = 0, Z = 0. If c x = we obtain the section which 
bisects the cone, so that instead of being hyperbolae the critical points fall on the lines 

1 r ( \V 2 " 

Y = m± l + (m-l)g 2 

m — 1 [ ^ ' 

Fig. 5 The projection of trajectories in the W = 0, A = subspace onto the X,Y plane, 
with V given by (2.12), for non-zero Cy (a) g 2 < m + 1; (b) g 2 = m + 1; (c) g^ 2 > m+1. 
The bold lines represent the same critical points as in Fig. 4. 

Fig. 6 The hemisphere at infinity for the 3-dimensional W = subspace with A 1 > 0: (a) 
gl 2 < 1; (b) gl 2 = 1; (c) 1< §1 2 < m + 1; (d) g^ > m + 1. 

Fig. 7 The hemisphere at infinity for the 3-dimensional W = subspace with A 1 < 0: (a) 
g 2 < m + 1; (b) g 2 > m + 1. 

Fig. 8 The plane Y = X or Y = (m+l)X/ (m— 1) in the Schwarzschild subspace. Both axes 
are 1-parameter families of critical points. The solution in the Y = X plane corresponds 
to the A = 0, A = limit of (2.20), while the solution in the Y = (m + 1)X/ (m — 1) plane 
corresponds to the limit g 1 = 0, Ai = of the solutions (B.13)-(B.21) in Appendix B. 

Fig. 9 The planes Y = m(X + k)/(m — 1) in the Schwarzschild subspace: (a) k = 0; (b) 
k < 0. The dashed trajectories are positive mass Schwarzschild solutions. If k > these 
trajectories are located in the X < —m\k\ region. 

Fig. 10 The hemisphere at infinity for the 3-dimensional Schwarzschild subspace: (a) with 
8, <f) defined by (2.24); (b) with 8, defined by (3.4). 
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